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Abstract 

Using the action principle we first review how linear density pertur- 
bations (sound waves) in an Eulerian fluid obey a relativistic equation: 
the d Alembert equation. This analogy between propagation of sound 
and that of a massless scalar field in a Lorentzian metric also applies to 
non-homogeneous flows. In these cases, sound waves effectively prop- 
agate in a curved four-dimensional "acoustic" metric whose properties 
are determined by the flow. Using this analogy, we consider regular 
flows which become supersonic, and show that the acoustic metric be- 
haves like that of a black hole. The analogy is so good that, when 
considering quantum mechanics, acoustic black holes should produce 
a thermal flux of Hawking phonons. 

We then focus on two interesting questions related to Hawking 
radiation which are not fully understood in the context of gravita- 
tional black holes due to the lack of a theory of quantum gravity. The 
first concerns the calculation of the modifications of Hawking radia- 
tion which are induced by dispersive effects at short distances, i.e., 
approaching the atomic scale when considering sound. We general- 
ize existing treatments and calculate the modifications caused by the 
propagation near the black hole horizon. The second question concerns 
backreaction effects. We return to the Eulerian action, compute sec- 
ond order effects, and show that the backreaction of sound waves on 
the fluid's flow can be expressed in terms of their stress-energy tensor. 
Using this result in the context of Hawking radiation, we compute the 
secular effect on the background flow. 



1 Introduction and summary 



According to General Relativity black holes are spacetime regions where the gravitational 
field is so strong that even light is unable to escape. It is indeed against all expectations- 
including his own- that Hawking discovered in 1974 [1, 2] that black holes are no longer 
black when Quantum Mechanics is taken into account. More precisely he showed that 
they thermally radiate (photons, neutrinos, gravitons, ...) with a temperature inversely 
proportional to their mass. This discovery had a major impact on theoretical physics. 

On a general level, Hawking radiation (HR in the following) finds its origin in the 
combination of the three pillars of physics as presently understood, namely Special Rel- 
ativity (with the speed of light c as fundamental constant), Quantum Field Theory (the 
relativistic version of Quantum Mechanics bringing h), and General Relativity (the metric 
and relativistic theory of gravity bringing Newton constant G). This triple origin can be 
seen in the expression for the temperature of a (uncharged and unrotating) black hole of 
mass M 

rp 1 

where kb is Boltzmann constant. 

An experimental detection of HR would therefore be of great importance. However, 
for a solar mass Black Hole (BH in the following), the temperature is only Th ~ 6x 10~ 8 K. 
This is much below the Cosmic Microwave Background temperature ~ 2.7 K. Hence there 
is no hope to detect HR from astronomical BHs which result from a gravitational collapse, 
since these should be more massive than the Chandrasekar mass (~ 3M sun ) [3]. 

Nevertheless, there is still hope to detect HR. First, light BHs could have been pro- 
duced in the early universe [4, 5]. If their mass is smaller than W 15 gr, their size would 
be of the order of one Fermi, and they will have enough time to explode in the present 
epoch at the end of the evaporation process caused by the HR. The associated burst of 
radiation would possess specific properties which could be detected in cosmic radiation 
[6, 7]. Unfortunately it is improbable to find such BHs if the universe underwent a period 
of inflation, as presently indicated by cosmological data (see for example [8]). Indeed, 
inflation would dilute them away like any other pre-existing population of objects. More- 
over there seems to be no efficient mechanism to produce them after inflation has stopped. 
Finally, at present, there is no observational evidence for excesses of radiation which could 
be interpreted as the bursts of radiation emitted during the last stage of the evaporation 
process. Nevertheless, the research in observational data is still currently active, see e.g. 
[9]. 

Perhaps a more promising experimental possibility of detecting HR, albeit indirect, 
arises from the astonishing analogy between light propagation in a curved spacetime and 
sound propagation in a (non-uniformly) moving fluid [10], or more generally, wave prop- 
agation in non-homogeneous condensed matter systems [11]. Indeed this analogy works 
so well that acoustic Black Holes (i.e. the region where sound is trapped) should emit 
a thermal flux of phonons for the same reasons that gravitational BHs emit HR. This 
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analogy is the central topic of this review article. Its precise content will be explained 
later. For the moment we further discuss the various ramifications of HR because this 
presentation illustrates the usefulness of the analogy. 

On one hand, HR solved the paradoxical aspects related to the proposal [12] made by 
Bekenstein in 1972 according to which black holes should possess an entropy proportional 
to the horizon area. Indeed the properties of Hawking Radiation confirm that black holes 
obey the laws of Thermodynamics and those of Statistical Mechanics (in fact not only BH 
possess thermodynamical properties but all event horizons in General Relativity, see e.g. 
[13]). In a sense, HR is a necessity dictated by the second principle of thermodynamics 
in the presence of horizons. We mention this aspect because it provides a first indication 
that HR should not rely on the details of the underlying fundamental theory, since its 
existence can be inferred by a thermodynamical argument. 

On the other hand, HR has raised, and it is still raising, new concepts and new ques- 
tions which concern intriguing (and poorly understood) aspects of space time dynamics, 
such as the microscopical origin of thermodynamics of horizons in the presence of gravity. 
In fact, the very existence of HR poses deep questions concerning Quantum Gravity and 
the validity of its semi-classical treatment based on quantum fields propagating in a clas- 
sical fixed background metric. It is this treatment which has been so far used to derive 
HR (see however [14, 15, 16]). Among the questions raised by HR, two shall be discussed 
in this article, in the light of the analogy with condensed matter systems. The first one is 
related to the occurrence of very high frequencies. This question is generally referred to 
as the trans-Planckian problem, and it has played an important role in the development 
of the above mentioned analogy [17, 18, 19]. The second issue concerns the backreaction 
effects, i.e. the consequences on the background metric induced by the non-linearities of 
the dynamical equations [20, 21]. 

1.1 Topics 

We first give the basis of the kinematical analogy between gravitational BH and sonic BH. 
Then we explain why an acoustic BH should also emit HR, and how the trans-Planckian 
question is addressed in the acoustic framework, knowing that the notion of acoustic 
wave ceases to exist when reaching the atomic scale. Finally we discuss the backreaction 
effects due to HR in the acoustic setting. Let us add that these last two sections should be 
considered as two first steps of a more ambitious program: to understand how to formulate 
backreaction effects at a deeper level when the microscopic structure of the fluid (or that 
of space-time) is taken into account. 

1.1.1 Near horizon propagation and Hawking Radiation 

It should be first stressed that both HR and the occurrence of very high frequencies result 
from the kinematical properties governing field propagation near a horizon. Therefore 



3 



these effects are not tied to BHs, nor even to General Relativity, but to near horizon 
propagation in general. When propagating backwards in time the field configurations 
giving rise to HR, there is a blue shift effect which grows without bound as the horizon is 
approached. Explicitly one finds 

fi(r) oc , (2) 

1 - r h /r 

where r — ru is the distance from the horizon, w the asymptotic frequency, and O(r) the 
frequency measured by freely falling observers crossing the horizon. This boundless growth 
results from the absence of a Ultra- Violet (UV) scale which could stop it. This absence 
is itself a direct consequence of Relativity: when studying fields obeying a relativistic 
field equation, their propagation in the near horizon geometry is governed by a simplified 
and universal equation which is scale invariant, the mass and the spin of the field having 
dropped out. This universal behavior guarantees that all fields are subject to HR, but 
also, commitantly, that no UV scale could possibly stop the blue-shifting effect of Eq. (2). 

This boundless growth is kinematical in the sense that it follows from the properties 
of outgoing null geodesies, which are the characteristics of the relativistic field equation. 
These characteristics hug the horizon, on both sides, according to 

x = r-r h = x e k{t - ta)/c , (3) 

where k is the so-called surface gravity of the horizon, and t is a time coordinate simply 
related to proper time at rest and far away from the hole. For a black hole of mass 
M, k = c 4 /4GM = c 2 /2rh (for an acoustic black hole, k is governed by the gradient 
of the velocity of the fluid evaluated at the horizon and t is the laboratory time). As 
we shall see in the text, this exponential behavior will be associated, both for acoustic 
and gravitational BHs, to the logarithmic singularity on the horizon of stationary modes. 
Moreover, in both cases, the regularity of the quantum state across the horizon implies 
the existence of HR with a temperature 

= (4) 

where c is the speed of light for gravitational BHs or the speed of sound for acoustic BHs. 
The important message is that there is a universal relationship between the characteristics 
obeying Eq. (3), the singular behavior of stationary modes subject to Eq. (2), and the 
existence of HR with a temperature given by Eq. (4) as a manifestation of regularity of 
the quantum state across the horizon. 



1.1.2 The trans-Planckian question 

The boundless growth of Eq. (2) is questionable [17, 22] both for acoustic and gravita- 
tional BHs. In the second case, this is because it has been obtained from quantum (test) 
fields propagating in a classical (passive) gravitational background. In other words, Eq. 
(2) follows from the implicit assumption that the quantum dynamics of gravity can be 
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neglected, an a priori good approximation for large BHs for which the curvature and the 
temperature are far from the Planck scale. 1 However, given the boundless growth, this 
approximation is questionable since the Planck scale is crossed in a logarithmically short 
time in the unit of (fc/c) _1 = 10~ 5 sec for a solar mass BH. 2 Therefore it is not unrea- 
sonable to assume that quantum gravitational effects will prevent the boundless growth 
of frequencies, thus invalidating the semi-classical treatment. However, if this aspect of 
the physics is not correctly described by this treatment, shouldn't we also question the 
existence and the properties of Hawking radiation, which have been derived using these 
approximate settings ? 

In the absence of a manageable theory of Quantum Gravity, no systematic analysis can 
be performed to settle this question (see [23, 25, 24] for attempts and for references). It is in 
this rather hopeless situation that the analogy with phonon propagation near an acoustic 
horizon is useful. The reason for this is rather simple. At low frequency, the analogy 
between photon propagation in a given gravitational background and phonon propagation 
in an acoustic metric is perfect: in both cases the quantum field obeys the d'Alembert 
equation. Therefore, if the acoustic metric possesses an event horizon, one expects to find 
a thermal emission of Hawking phonons. However, when reaching the molecular scale d c 
(a short distance with respect to the characteristic length scale = (fc/c 2 )^ 1 ), phonon 
propagation is modified and no longer described by the d'Alembert equation. Indeed, in 
fluids, there is a UV regime (for wavelenghts smaller than d c ), in which the phonons are 
described by a non-linear dispersion relation, which breaks the Lorentz symmetry which 
emerges at low frequency. 

1.1.3 Two lessons from Acoustic BHs 

This breakdown of Lorentz symmetry raises the interesting question of whether HR will 
be modified, given that Hawking phonons emerge from configurations with wavelength 
reaching the molecular scale. The outcome of the analysis is that HR is robust [18, 19]. 
That is, to leading order in the parameter K,d c >C 1, the properties of HR, thermality and 
stationarity, are preserved. From this we learn the first lesson: the properties of HR are 
not tied to the relativistic (scale invariant) mode equation. 

Besides this robustness, there is another important lesson which is delivered by the 
analogy. One notices that the unbounded growth of frequencies associated with Eq. (2) 
is exceptional, in that it is found only for the relativistic dispersion relation uj 2 = c 2 p 2 . 

1 The Planck curvature is defined as lp 2 , where lp = \J Gfi/c 3 = 1.616 x 1CP 33 cm, and 
the Planck temperature as E p /k b , where E P = sJhcFjG = 1.221 x 10 19 GeV\ 

2 One might question this statement since one can think that one can avoid large fre- 
quencies by a judicious choice of coordinates. However this is not the case: the gravita- 
tional tree-level scattering amplitude between an outgoing wave and an infalling one is a 
scalar, and grows boundlessly when approaching the horizon as fi( r ) does it in Eq. (2) 
[23, 24]. 
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For all non-linear dispersion relations, the new UV scale p c = l/d c effectively cuts off the 
growth of frequencies. This generic behavior raises doubts about the validity, not on the 
existence of HR, but on the free field dispersion relation lo 2 = p 2 , used to characterize field 
propagation in the near horizon region [26], see also [27] for a broad audience presentation 
of this point. Indeed the use of this relation amounts to assume that no UV scale would 
modify the two-point correlation function when taking into account short distance effects. 

1.1.4 Backreaction effects on acoustic metrics 

It should be emphasized that HR has been derived from a linearized treatment, in which 
the curved metric entered only as a passive background. However, in a full treatment, 
the degrees of freedom of the metric (or those at the origin of the metric) should also be 
treated dynamically and quantum mechanically. The first step to include the backreaction 
is clear and rather standard. It corresponds to the usual secular effects, which occur to 
second order in perturbation theory. These secular effects appear classically and quantum 
mechanically, both for gravitational BHs and for the acoustic BHs. It is therefore of 
value to compare how these influence both kinds of systems since, a priori, the analogy 
works only at the kinematical level. Indeed, the evolution of a gravitational black hole is 
determined by Einstein equations, whereas the dynamics of a sonic black hole is governed 
by hydrodynamical equations. 

1.1.5 Plan of the paper 

The paper is organized as follows. 

In section 2 we describe the classical propagation of sound waves in inhomogeneous fluid 
flows. This provides the formal structure which governs the hydrodynamical-gravitational 
analogy. 

In section 3 we recall the basics features of quantum field theory in a curved metric. 
In section 4 we show how HR emerges when an acoustic BH forms. 

In section 5 the trans-Planckian question is addressed when taking into account the mod- 
ification of the dispersion relation at short distance. 

Section 6 is devoted to the study of the backreaction on the fluid flow due to the energy 
momentum tensor of the quantum phonon field. 
In section 7 we conclude with some remarks. 

To further appreciate the analogies and the differences between the gravitational and hy- 
drodynamical cases, we have added for non-expert readers an extended Appendix wherein 
we review the major features of BHs in General Relativity. 
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2 Sound waves in inhomogeneous flow 



As discussed in the Introduction, Hawking Radiation is not related to the dynamics of 
the curved metric. One may therefore expect to find HR in different physical contexts 
whenever the propagation of a quantum field is similar to that in a BH metric. The simplest 
example is provided by a fluid whose motion becomes supersonic in some region. In a 
seminal work, dated 1981 [10], Unruh showed that sound waves in this system propagate 
exactly as light does in a black hole metric. For reviews on this topic we refer to [28, 29, 30]. 
In this section we analyze how this result is obtained from a linearized treatment. We 
then study the next order effects (backreaction). 



2.1 The action principle 

The action principle constitutes the most straightforward way to derive the equations of 
motion of a system (Euler-Lagrange equations) and the conservation laws associated to 
its symmetries (by making use of Noether theorem). This procedure can also be applied 
to fluid mechanics. 

Let us consider a fluid which is irrotational (i.e. the velocity v can be expressed as gradient 
of a potential ip) and homentropic (i.e. the pressure P is a function of the density p only). 
In the absence of external forces, the Eulerian equations of motion for this fluid can be 
derived from the action [31, 32, 33] 



pip + ip(Vi/)) 2 + u(p) 



(5) 



where u(p) is the internal energy density and a dot means derivative with respect to 
newtonian time t. 

Variation of S with respect to ip gives the continuity equation 

p + V • (pv) = , (6) 
with v = Vip. Variation of S with respect to p yields Bernoulli equation 

ip + ^v 2 + ^(p) =0, (7) 
with p(p) = The gradient of the Bernoulli equation gives the Euler equations 



ir+ (v- v)v+ ivp = o , (8) 



where the pressure P = J pdp. 

The invariance of the action under translations 



ip{x l ) — > il>(x % — a % ) , 

p{x i ) - p{x-a*), (9) 
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leads, through Noether theorem, to the momentum conservation law 

d t {pVi) + V • (vvi) + pd % p = , 

or, equivalently, 



(10) 



dt(pvi) + djUij =0 , (11) 

where Tlij = pviVj + &ijP is the momentum flux tensor and pVi the momentum density. 
Finally, the invariance of the action under rescaling of the velocity potential, tp — » if) + a, 
brings again the continuity equation (6). 



2.2 The Acoustic Metric 

To obtain the propagation of sound waves, we expand the action S up to quadratic order in 
the linear fluctuations pi and ip\ around some background configuration po and Vo- This 
configuration describes the mean flow, namely a solution of the (unperturbed) equations 
of motion (6, 7). Therefore, the linear terms in pi and ?/>i drop out, and one has 



S — So + S2 , 



where 



S 2 



J- 



c 2 1 

piV'i + - — pi + ttMVV'i) 2 + Piv ■ VV>i 

IpO 2. 



(12) 



(13) 



In the latter equation, we have introduced the speed of sound c, which is defined by 
c 2 _ p dji p Qr g j m pij c jty 5 m this paper we shall consider c = const. 

One can eliminate pi in favor of tpi by using the pi equation of motion 

„2 



ipi + v ■ Vl -\ pi 

po 



0, 



and find 



ipotWif-g^i+iT-VVi) 5 



(14) 



(15) 



This action gives the evolution equation for the perturbation of the velocity potential ipi : 
-d t [^(dt^+vo- V^i)] + V-[v [-g(<9 t Vi+tfo-VVi)] +PoV^i} = . (16) 

Eqs. (14) and (16) describe the propagation of linear fluctuations, i.e. sound waves. 

The crucial point for us is that Eq. (16) can be rewritten in a four dimensional 
notation as 

d M (,rU^i)=0, (17) 



where 



Moreover, if we set 



-1 -vl 



-99^ 



(18) 
(19) 
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then Eq. (17) corresponds to the d'Alembert equation in the curved metric described by 
g"" [10, 34, 35], namely 

-jq-d^J—gg^d^x) = D^i = . (20) 

Similarly, the action Si of Eq. (15) can be rewritten as 

Si = ~\ J d^x^—gg^d^d^x . (21) 

One can also introduce the covariant "acoustic metric" which is the inverse of 
and reads 

/ i 1 2\ -T 

Po -(c -v ) -v 



9** = V . " r - (22) 

c y -«o 1 

where 1 is the identity 3x3 matrix. The corresponding "acoustic line element" is 

ds 2 = g^dx^dx v = ^^-{c 2 -Vo)dt 2 ~28 ij vl)dx : >dt + Si j dx i dx : >^ . (23) 

In brief, according to Eq. (20), sound waves propagation coincides with the propaga- 
tion of a minimally coupled massless scalar field in the curved spacetime described by the 
acoustic metric g^v? Unlike the particles constituting the fluid, which see spacetime as 
Galilean and flat, 4 sound waves propagate in a four dimensional curved metric <? M „ (the 
origin of the curvature is rooted to the non homogeneity of the flow, see [11]). As we shall 
see in the following subsection, this implies that sound propagates along the null geodesies 
of g^u exactly as light moves along the null geodesies of spacetime. This analogy between 
propagation of sound waves and propagation of massless radiation in a gravitational field 
allows to translate, using the language of differential geometry, geometrical features of 
General Relativity into the hydrodynamical model. 

It should be also noted that the three dimensional spatial sections t = const, of the 
acoustic metric g^v are conformally flat. When considering an inward radial flow, the 
acoustic metric <? M „ resembles to that of a Schwarzschild metric expressed in the so-called 
Painleve-Giillstrand (PG) coordinates, for which the spatial sections are also flat (see Eq. 
(199) in the Appendix C). Note that the resemblance is up to an overall conformal factor 



3 The general action for a non-interacting scalar field 4> propagating in a curved space- 
time has the form [36] 



S = - I d*Xy 



(24) 



where m is the mass of the field, and £ measures the non-minimal coupling to the Ricci 
scalar. Eq. (21) corresponds to the minimal coupling massless case, namely m = £ = 0. 

4 Note that in [37] Moncrief started with a four dimensional relativistic space (flat or 
curved), considered a relativistic flow and derived the equations for the linear perturba- 
tions. In this case as well, one finds that the equations of motion are governed by a four 
dimensional "acoustic metric". 
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which, as well known from differential geometry, leaves unchanged the equations yielding 
the null geodesies. Hence, in our case, the conformal factor does not affect the propagation 
of sound waves [3, 35]. 



2.3 Geodesies from geometric acoustics 

In this section, we translate in the acoustic setting the derivation of the geodesic equations 
of General Relativity given in [3] and in Appendix A. 

Let A be the typical wavelength of sound waves and L the typical length over which the 
amplitude and the wavelength of the waves vary. Geometric acoustic is valid whenever 
A< L. 

Focussing on waves that are highly monochromatic in a region ~ L, we split the velocity 
potential t/>i into a rapidly changing real phase S and a slowing varying amplitude a 

tpi = 5R{ae lS } . (25) 

Introducing the wave vector = d^S, sound rays are curves normal to the constant 
phase surface S = const. The differential equations for sound rays is dx^ /ds = where 
s is a parameter. Using Eq. (20), one gets 

g^K»K v = , (26) 

i.e. if M is a null vector. After covariant differentiation of the equation above, one can 
show that K 1 * satisfies the geodesic equation (see Eq. (174) in Appendix A). Thus, sound 
propagates along the null geodesies of the four dimensional acoustic metric (22). 

It is interesting to re-express the above equation in a 3 dimensional language. Ex- 
panding the null condition, and using the Newtonian time in place of s, we get 

g lu ,K»K" = g lu ,*g*£ = 

= § [~(c 2 - ~ 2*««5^ + = o • (27) 

This equation implies that 



dx 

dt- V ° 



c, (28) 



where the norm is taken in Euclidean space. This is the usual (non-relativistic) equation 
describing sound propagating with velocity c with respect to the medium moving with 
velocity vo- 



2.4 Pseudo energy-momentum tensor 

The energy-momentum tensor is best defined as the functional derivative of the relevant 
action with respect to the background metric (see for example [38]). In the case at hand, 
two metrics are at our disposal: the real Minkowski spacetime metric and the acoustic 
metric defined in Eq. (22). Differentiation with respect to the former leads to the real 
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energy and momentum (see Eq. (11)), while differentiation with respect to the acoustic 
metric leads to the so called pseudo energy-momentum tensor (for a detailed discussion 
see [33]). 

Writing the action £2 in terms of the acoustic metric as in Eq. (21), the pseudo 
energy-momentum tensor for the fluctuations is 

2 SS 1 

Tu„ = 1=-; = duipid v ipi - -g^d a ipid a ipi ■ (29) 

V — 9 "9 2 

Its covariant components read 

T ti = fadah.--^ [(i/>i + w>i-i7 ) 2 -c 2 (vvi) 2 ] , 

Tij = di^pid^! + ^Sij [(V>i+VV>i-iTo) 2 -c 2 (VV>i) 2 ] • (30) 

Because of the invariance of S2 under diffeomorphisms, this tensor satisfies the "covariant" 
conservation law 

-^LcV {V=gT» a ) + r^r^ = , (31) 

where the connection coefficients T^ v are defined as 

= \g afi (SW,* + 3^/3,m " 9^,fj) ■ (32) 

To appreciate the difference between the true energy- momentum tensor and T M ", as de- 
fined by Eqs. (30), and to understand the meaning of Eq. (31), one should note that for 
our system there are two distinct operations which can be called "translation" in the (say) 
^-direction: 

a. the fluid, together with any disturbance on it, is translated in the ^-direction; 

b. the fluid itself is left fixed, while the disturbances are translated in the 2-direction. 

The first operation is a symmetry of the action provided the background space is ho- 
mogeneous. The associated conserved quantity is the true newtonian momentum J d 3 xpv. 
Operation b instead is a symmetry only when both the background space and the fluid 
are homogeneous. In this case the T™ v vanish since the acoustic metric is constant and 
Eq. (31) yields the conservation of pseudo energy-momentum. 

In the general case, r™„ 7^ and Eqs. (31) describe how in an inhomogeneous flow field 
energy and momentum are exchanged between the waves and the mean flow [33]. 

2.5 Acoustic black holes 

As we have seen in subsection 2.3, sound propagates along the null geodesies of the acoustic 
metric g^v of Eq. (22). The latter shows deep resemblances with the Painleve-Giillstrand 
form of the Schwarzschild metric (see Appendix C Eq. (199)). 

For simplicity, we shall not consider radial flows with spherical symmetry. We consider 
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instead linear flows in the z direction, and assume no dependence on the perpendicular 
directions x and y. We shall also restrict our attention to the null geodesies which are 
parallel to the flow, i.e. dx = dy — 0. According to Eq. (27), those obey 

= [-(c 2 - v 2 )dt 2 - 2vdzdt + dz 2 ] . (33) 

When supposing that v < 0, i.e. the fluid moves from right to left, the sound waves 
propagating downstream follow 

|=,-c, (34) 
whereas the sound waves propagating upstream obey 

f t =v + c. (35) 

Now, if the fluid motion becomes supersonic, the surface \v\ = c separates the up- 
stream null geodesies into two classes. In the region where the motion is subsonic \v\ < c, 
upstream waves can propagate to the right (^| > in Eq. (35)). Instead, in the region 
where the motion is supersonic, upstream propagating waves are dragged by the fluid mo- 
tion to the left (^| < in Eq. (35)). The surface of the fluid where \v\ — c represents the 
sonic horizon. It can be crossed only by downstream propagating sound waves. Inside this 
surface, i.e. in the supersonic region, sound cannot escape: it is trapped by the supersonic 
motion of the fluid in the same way as light is trapped inside the horizon of a black hole 
by the gravitational field. For these reasons, it is appropriate to call this configuration an 
"acoustic black hole" or a "dumb hole" [10] (dumb = unable to speak). 
For stationary flows one can introduce null coordinates x ± 



Thus, upstream (outgoing) waves follow x~ = const., whereas downstream ones (ingoing) 
follow x + — const. 

A diagram representing the propagation of these waves is given in Fig. 1. One should 
appreciate the strict resemblance between this diagram and the one depicting the behavior 
of light rays in the Schwarzschild black hole geometry (see Fig. 5 in Appendix C). 
It is interesting to further analyze the behavior of the geodesies in the near horizon region, 
because it determines the properties of Hawking radiation when studying quantum effects. 
In the near horizon region one can develop the (stationary) velocity profile as 

v{z) = -c + ck(z- z h ) + 0{z- z h f , (37) 

where zk is the location of the sonic horizon and we have introduced k which is given by 
k — k/c 2 , where k is the "surface gravity" of the sonic horizon, defined as [10, 34] 



k = - 



1 d(c 2 -v 2 ) 



2 dz 



(38) 
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z 



Figure 1: Propagation of sound waves in an acoustic BH. Families of in- 
falling (downstream) x + = const, and outgoing (upstream) x~ = const, null 
geodesies are depicted. The verticallgxis is the laboratory time t, and the 
horizontal one is the z axis. The fluid is accelerating and flowing towards 
the left. It becomes supersonic at z = Zh which is represented by the thick 
vertical line, the horizon. On its left upstream propagating waves are swept 
away by the supersonic flow and propagate inwards. 



Inserting the near horizon expression of v(z) in Eq. (35), one sees that outgoing 
geodesies (x~ = const.) leave the horizon following an exponential law governed by k 

z(t) -z h = Ce CKt . (39) 

This implies that the frequency (measured in the rest frame of the atoms) of an outgoing 
sound wave emitted near the sonic horizon will suffer an exponential Doppler shift 

Q(t) ~ tt e- CKt , (40) 

as it escapes from the hole. (The proof will be given in subsection 4.2.2, see Eq. (99)). 
Comparing these expressions with those of a radiation field in a gravitational black hole 
metric (Eqs. (195, 197) in Appendix C) one verifies that there is a complete analogy 
between light propagation in black hole geometry and sound propagation in an acoustic 
black hole geometry. This should no longer be a surprise since the field actions and the 
metrics are identical. 



2.6 De Laval nozzle 

Gravitational black holes are formed by the collapse of very massive objects (M ~ 3M aun ). 
Sonic black holes on the other hand do not require such extreme conditions for their 
realization. They can be formed by manipulating laboratory size systems. 
A simple realization is provided by a converging-diverging nozzle, called a de Laval nozzle, 
well known in jets engine engineering (see for ex. [39]). It is schematically depicted in 
Fig. 2. The fluid flows from right to left and the nozzle is pointing along the z-axis. Let 
the cross section area be A{z). We assume that the transverse velocities (in the x and y 
direction) are small with respect to the velocity along the nozzle axis (i.e. z). We have 
then a quasi one dimensional flow which we further assume to be stationary. 
In this case the continuity equation yields 

p(z)A(z)v(z) = D , (41) 

where D is a constant. The Euler equation (8) reads 

dv c 2 dp 

v— + --f=0. 42 

dz p dz 

By differentiating Eq. (41) to eliminate p in Eq. (42), we obtain 

where a prime ' means derivative with respect to z. For subsonic flows (i.e. |v| < c) the 
velocity decreases as the area increases and vice versa as everyday experience tells us. At 
the waist of the nozzle (say z = 0), A' = and therefore v' = 0. 

Suppose the fluid starts from z = +oo with some velocity, say u+oc. The velocity increases 
through the converging nozzle (i.e. A' /A < 0, v' /v > 0) reaching the maximum value 
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Figure 2: A de Laval nozzle. As explained in the text, it is possible to set 
the system to have a regular flow that is subsonic on one side of the waist 
and supersonic on the other. The acoustic metric associated with this flow 
is that of a one-dimensional black hole whose null geodesies are represented 
in Figure 1. The sonic horizon (\v\ = c) is located at the waist of the nozzle. 

| «o| < c at the waist, and then decreasing to the value «_oo along the diverging part of 
the nozzle (A' /A > 0, v' /v < 0). 

The picture changes completely when the motion is adjusted to become supersonic. Prom 
Eq. (43) one sees that the transition from subsonic to supersonic flow can occur only 
at the waist of the nozzle (unless the acceleration vv' diverges when |v| = c, which we 
consider as unphysical). Flowing along the divergent part of the nozzle (z < 0) the fluid 
increases its velocity since, for |v| > c, A' and v' have the same sign in Eq. (43). This 
non intuitive behavior is a well known feature of supersonic fluid motion. The increase 
of velocity is associated with a decrease in the density (and pressure) in the supersonic 
region (z < 0). 

The above fluid configuration therefore corresponds to that of a sonic black hole with the 
sonic horizon located at the waist of the nozzle 2 = 0. The trapped region is the supersonic 
region (|v| > c) on the left of the waist. 

When the velocity of sound is constant, the nozzle equation (43) can be immediately 
integrated yielding for the acoustic black hole configuration 

A=^A H e^-^ 2 , (44) 
|v| 

where Ah is the area of the waist of the nozzle, and where integration constants are chosen 
such that for A = Ah we have v = — c. It is interesting to compare this with the subsonic 
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solution (\v\ always < c). In that case, one has 

^^He^. (45) 



In both cases, the Bernoulli equation gives 



^+c 2 ln-^=0, (46) 
2 po 



p = poe-* 2 / 2c2 , (47) 

where the relation p = <? lnp/po, following from the constancy of c, has been used. 

The acceleration of the particles of the fluid is given by a — vv' . In terms of this parameter 

the nozzle equation (43) can be rewritten as [34] 

v 2 c 2 A' 
a ~ c 2 —v 2 A ' 

and on the horizon (using de l'Hopital rule) we have 



(48) 



2 

Cijf = C 



A><_ 
~A 



(49) 

H 

This definition coincides in our case (c = const.) with the surface gravity defined in Eq. 
(38). 

We have seen that the fine tuning condition A' = forces the acoustic horizon to form 
and remain exactly at the waist of the nozzle. The flow will in fact self-adjust to satisfy 
this fine tuning which keeps the acceleration finite at the acoustic horizon. On the other 
hand, in presence of an external force F entering Euler equations (8), the nozzle equation 
(43) gets modified, namely 

A v \ c 2 ) pc 2 

The horizon location now satisfies the equation A'/A\h = — F/{pc 2 )\n and does not 
coincide in general with the waist of the nozzle. 

Finally, the formation of a sonic black hole can be obtained by a continuous deforma- 
tion of the velocity profile. Starting from an initial subsonic configuration, described in 
Eq. (45), one obtains the black hole configuration described in Eq. (44) by lowering the 
pressure in the exhaust region (large negative z). This evolution is depicted in Fig. 3 (see 
[40]). 



2.7 Backreaction of sound waves 

Because of the non-linearity of the hydrodynamical equations, linear perturbations will 
backreact on the mean flow. The dominant effect is the secular change that they induce. 
The backreaction of the waves on the underlying fluid can be evaluated following the same 
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Figure 3: Spacetime diagram representing the formation of a sonic black 
hole. The solid line denotes the sonic horizon. At late times, the acoustic 
metric coincides with the static one of Fig. 1 as indicated by the trajectories 
followed by outgoing null geodesies. 
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procedure used to calculate the backreaction due to the gravitational waves in General 
Relativity, as it can be seen from Appendix B and [3]. One writes the physical quantities 
which describe the fluid as a term describing the mean flow plus a fluctuation, i.e. 

p = PB + pi ■ (51) 

The subscript B stands for backreaction, i.e. these quantities include the backreaction of 
the waves. The functions pi and ipi are linear in e, where e is a dimensionless parameter 
which characterizes the amplitude of the perturbations. The continuity and Bernoulli 
equations are then expanded to order e 2 

p B +pi + V- [{p B +pi)(v B + vi)] = , (52) 



; ; 1 ,_ _ . 2 , n du 

Vb +Vi + -{VB +vi) + p(pb) + 



PB 



i d' 2 p 
P1 + 2W 



pi = . (53) 



Since pB and vb contain terms of order e 2 only, the terms linear in e must vanish, thus: 

pi + V • (pbvi + piv B ) = , (54) 
c 2 

V>i + vb ■ vi H pi = . (55) 

pB 

As we have seen in subsection 2.2, these lead to the wave equation for a massless scalar 
field ipi. The novelty is that now the d'Alembert operator O g is evaluated in the "mean" 
acoustic metric determined by pB and vb- 
Prom Eq. (22), we have 

„b _ Pb ( -(c 2 - v%) -v% 



C ^ -VB 1 J 

Taking into account Eqs. (54, 55), the mean flow (pb,vb) satisfies the backreaction 
equations Eqs. (52, 53) without the linear terms in e: 

PB +V- (p B v B ) + V- (pivi) =0 , (57) 

1 1 c 2 

V>s + 77«s + M(Ps) + t:v( - - — p{ = . (58) 

2. I IpB 

There is an elegant way to rewrite these equations in terms of the four dimensional 
pseudo energy-momentum tensor of the fluctuations (see Eq. (29)), now evaluated with 
respect to the mean acoustic metric g^ v . We first note that 

i PB I '. _* \ i 

pivi = 7T (ipi + v B ■ «i)«i 

= -^f( T « + "b^'j) = V^T? , (59) 

and 

1 ,2 C 2 1 PB „ , . 

2 Wl -2^ Pl = -2- T ' (60) 
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where T = g% T a g is the trace of the pseudo energy momentum tensor. Then Eqs. (57, 
58) become 

PB + V-(pBV B )+d 1 (\/ r ^T°) = , (61) 

i>B + \tl+n{p B )-\?fT = Q . (62) 

Taking the gradient of Eq. (62) one obtains the Euler equation which determines the 
dynamics of the mean flow 

i>B+VB-VvB + -VP B -lv(—T)=0. (63) 
p 2 \ c I 

Notice that if c 7^ const, extra terms will appear in Eq. (62, 63). 

There is yet another concept which should be introduced. One often deals with an en- 
semble (stochastic or quantum) of fluctuations rather than with a single superposition of 
fluctuations. Then, to write the backreaction equations for the mean flow, one should 
replace T M „ in the above equations by its expectation value (T^) in the ensemble under 
consideration. These expectation values are functions of Pb,vb through Eqs. (54, 55), 
therefore Eqs. (61, 62) are indeed the dynamical equations governing the evolution of the 
mean fields Pb,vb, in the Hartree-Fock approximation. 

In conclusion, one sees that the continuity equation receives a contribution from the 
pseudo momentum density, while the Bernoulli equation gets a contribution to p propor- 
tional to the four dimensional trace of the pseudo energy momentum tensor. 



3 Quantum Field Theory in curved spacetime 

As we have seen the propagation of sound waves in a non homogeneous fluid flow corre- 
sponds to that of a massless scalar field in a curved spacetime described by the acoustic 
metric. Hence quantization of these sound waves should be discussed within the frame- 
work of QFT in curved spacetime (see e.g. [36]). We outline here this procedure. 
We start by recalling the canonical quantization of a scalar field <f> in Minkowski spacetime 
in the Heisenberg representation, where the metric tensor is 7/ M „ in cartesian coordinates. 
The action of <j> is 

S = ~\J ^xn^d^d^. (64) 
Hence (f> satisfies the massless Klein-Gordon (KG) equation 

□0 = rf'd^ = . (65) 

To quantize <j>, it is appropriate to separate the positive and negative frequency solutions 
with respect to the Minkowski time t. That is, one writes the operator 4> as 

<j>{t, x) = [ojWj (t, x) + a\u* (t,x)], (66) 
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where 

— Uj(t,x) — ~iujjUj{t,x) , ojj > . (67) 

The positive frequency solutions Ui(t,x) and their complex conjugates form a complete 
orthonormal basis with respect to the Klein-Gordon inner product defined as 

U1J2) = -ij d A x(hd t ti - fidth) , (68) 

where /1 and /2 are solution of the KG equation. The Fock space is constructed from the 
vacuum state |0) annihilated by the destruction operators a; 

a,i\0) = . (69) 

The one-particle states are obtained by acting with the creation operators a] on this 
vacuum. Similarly, the many-particle states are constructed by acting repeatedly with the 
creation operators. 

Choosing as the orthonormal basis the plane wave modes 

_ 1 e -iw k t+ik.s ^ (7Q) 

Vl67T 3 U!fc 

where Wk = \k\, the commutation relations for the creation and annihilation operators are 

[a % ,ai kl \ = h5\k-k') , (71) 

[a % ,a p ] = Q=[al,a\] . (72) 

This procedure can be extended to curved spacetimes described by a metric tensor 
g uv - The action for a minimally coupled massless scalar field propagating in a curved 
spacetime is given by Eq. (24) and the field equation reads O g <f> = 0, where D g is the 
d'Alembert operator in the metric g uv - 

In the canonical quantization procedure, a crucial difference with respect to flat spacetime 
is that in curved spacetime there is no global definition of time and, therefore, there is 
no natural splitting of modes between positive and negative frequency solutions. Differ- 
ent choices of positive frequency solutions lead, in general, to different definitions of the 
vacuum state [36]. 

When there is a timelike vector field £ M leaving invariant the spacetime metric, one recovers 
a notion of positive frequency solutions: 

£ M V M Uj = —iujjUj , Wj > . (73) 

When the spacetime is not stationary, one looses the above criterium to define a single 
set of positive frequency modes. However, when the spacetime possesses asymptotic sta- 
tionary regions in the past and in the future, two sets of positive frequency modes can be 
defined. Let us call u] n the orthonormal set of positive frequency modes in the remote 
past. We can then expand the field <j> as 

1 \ y r in in , ini in*i i n A\ 

<t> = u * + a i u i 1 ■ ( 74 ) 
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Similarly, let u° ut denote the set of positive frequency modes in the late future. The 
corresponding expansion for our field is then 

i \ v r out out . out] out*-\ 

4> = 2__ l [a i Ui 'Uj ] . (75) 

i 

Each set of modes defines its own vacuum state: \in) and \out) respectively. They are 
defined by 

at\in) = 0, Vi , (76) 

and 

a° ut \out) = 0, Vi . (77) 
Considering our quantum field to be in the in vacuum state \in) (i.e. a state with 
no particles in the past), the point is to investigate the physical effects induced by the 
time dependence of the background metric at late times, when the geometry has settled 
down again to a stationary configuration. In general, the state \in) will be not perceived 
as a vacuum state by particle detectors situated at late time. These detectors will indeed 
react to the presence of out quanta. To quantify these effects it is sufficient to compute 
the Bogoliubov transformation which relates the two sets of modes [36]. 

Since both sets of modes are complete, one can always expand modes (and operators) 
of one set in terms of those of the other set: 

out \ / in . 

u i = l^WiU-i + (3jiUi ) , 

i 

out \ ^ / * in r>* in t \ / -rn\ 

i 

These relations are the Bogoliubov transformations and the matrices o^-, , /3ij are called 
Bogoliubov coefficients. They are defined by the overlaps 

/ out in\ rt / out in* \ /i-?n\ 

a ij = {u i ,Uj ) , fiij = -{Ui ,u 3 ), (79) 
where ( , ) denotes the Klein-Gordon product extended to curved spacetime 

(/i,/ 2 ) = -i J dsruidpf; - fidpfi) . (so) 

The symbol E designates a Cauchy hypersurface for the d'Alembert operator. The dif- 
ferential dS M = dEn M , where d£ is the volume element and n M is a future directed unit 
normal vector to E. 

The non- vanishing character of the 0ij coefficients in Eq. (78) means that the positive 
frequency modes u™ are not positive frequency in the out region. This implies that the 
vacuum state \in) contains out quanta. Indeed, using Eq. (78), the mean occupation 
number of the i th out quanta, N° ut = a° ut] of', is 

(in\Nr\in)=J2\^\ 2 • ( 81 ) 

3 

Hawking radiation is an application of this result. In that case, the difference between the 
\in) and \out) vacua is a consequence of the formation of the black hole. 
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4 Hawking radiation in acoustic black holes 



4.1 Introduction 

Given the equivalence between the action governing massless radiation propagating in a 
black hole geometry and that of sound waves in an acoustic black hole geometry, one 
expects, when sound waves are quantized, to obtain Hawking radiation in the form of a 
thermal emission of phonons when a sonic hole forms [10]. We shall outline here how this 
result is obtained, giving a more detailed derivation in the next section. 
The phonon field ipi is quantized exactly like the massless scalar field <j> in the former 
section. Its dynamics is governed by the action of Eq. (21) and it obeys the usual 
commutation relations. 

Let the acoustic metric be the one associated to the flow (from right to left) of 
a fluid in a de Laval nozzle, as discussed in section 2.5. Suppose the initial configuration 
(t — > — oo) corresponds to a stationary subsonic flow. In this case, the normal modes 
of the field are asymptotically plane waves e~ lClx e~ lVtx / c , where x ± are defined in 
Eqs. (36). The former modes propagate downstream, from right to left, and are Doppler 
shifted to shorter wavelengths when the velocity of the flow increases. The latter propagate 
upstream, from left to right, and are Doppler shifted to larger wavelengths in the same 
situation. Both modes are positive frequency with respect to the Newtonian time t, and 
stay positive frequency as long as the flow stays stationary and subsonic. 
Expanding the field in this basis, one defines the in vacuum \in) which contains no quanta 
(i.e. no phonons) at t = — oo. 

Now suppose that the velocity profile varies till an acoustic horizon develops, and 
then settles so as to engender a stationary black hole geometry of the form described by 
Eq. (22). When the horizon forms, the modes propagating upstream separate into two 
classes, as clearly seen in Figure 3. Those on the left of the horizon remain trapped in 
the supersonic region. Instead those on the right succeed in escaping the hole and suffer 
a Doppler shift which becomes infinite, as one might guess from Eq. (40). Just before 
the formation of the horizon, one can show that the Doppler shift exponentially increases. 
This means that an initial plane wave of frequency SI behaves at late time and far away 
from the horizon as e~ lSlx , where X~ scales as 

X~ oc -e~ KX ~ (82) 

with k — <?k the surface gravity of the sonic horizon. 

These distorted modes are no longer positive frequency with respect to the time t. In other 
words there is a non-trivial (f3 =fc 0) Bogoliubov transformation between these in modes 
and the stationary out modes ~ e~ %ut of positive frequency to. When evaluating these 
Bogoliubov coefficients, one obtains, following Hawking [2], that the mean occupation 
density (i.e. per unit frequency and per unit time) of outgoing asymptotic quanta obeys 
the Planck distribution: 

Wlm) = e ^ /K 1 BTH _ 1 . (83) 
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where 

kbT„ = ^- . (84) 

An estimate for the order of magnitude of this temperature gives (see Visser's contribution 
in [30]) 

T n 1 r 1 rl(r — ,A~l 

(85) 



1.2 x lO^Kmm 



1 d(c-v) 
c dz 



1km/ sec 

For the case of a supersonic flow of water through a 1mm waist we have T H ~ 10 K, 
while for liquid Helium T H ~ 10~ 4 K. In the case of Bose-Einstein condensates this 
temperature is ~ 10 nK, which is not so low if compared with the usual condensate 
temperature ~ 100 nK (see [41]). 

4.2 Formal derivation 
4.2.1 Introduction 

We give a detailed derivation of the Hawking Radiation in acoustic black holes in the usual 
relativistic way [42, 43] in a manner which on one hand reveals the essential ingredients 
of this radiation, and on the other hand will show why the properties of HR are robust 
when using a nontrivial dispersion relation which breaks Lorentz invariance. 
In order to keep the discussion as simple as possible, we shall adopt several simplifications 
which do not affect the key features of HR. 

As in section 2.6, we assume the acoustic black hole is described by a quasi one- 
dimensional flow of the fluid along the axis of the nozzle (the z axis). The corresponding 
acoustic metric is 

ds 2 = ^ [~c 2 dt 2 + (dz - vdtf + dx 2 + dy 2 ] , (86) 

where the relevant variables, v = v z and p, depend only on t and z. 

In this metric, the phonon field if) is described by the action, see Eq. (21), 



(87) 



HR arises from the near horizon propagation of the modes of the quantum field V- In this 
region (z — z^ -C 1/ ' k) p varies smoothly, so we approximate it as a constant in the sequel. 
This approximation implies neglecting of the potential barrier associated to a varying p 
which causes an elastic scattering of the modes. This would reduce the asymptotic flux, 
since the transmission coefficient across the barrier is smaller than unity, without affecting 
the very essence of HR. From now on, we shall work with a rescaled field 4> = */p~ip- 

For further simplicity, we shall consider only waves that do not depend on the trans- 
verse variables x, y. (A transverse momentum can be easily included as it will be conserved 
and will act as a mass square in the equation for the longitudinal part of the modes). In 
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brief we shall study the solution of the reduced two-dimensional equation obtained from 
Eq. (87), namely 

[{d t + d z v(t, z)){dt + v(t, z)d z ) - c 2 d 2 ] <f>(t, z) = . (88) 

In the next subsection we study the stationary case wherein v does not depend on time. 
This is in anticipation to the fact that Hawking Radiation is a stationary process. 

4.2.2 The stationary modes 

Because of time translational invariance, it is appropriate to decompose the scalar field as 

<t>(t,z) = J duj ^(z)e- iut + h.c^j . (89) 

The stationary modes <t>w{z) obey, see Eq. (88), 

(w + id z v)(u + ivd z )<f> w = -c 2 d 2 <f> w . (90) 
The solutions are orthonormalized by the Klein-Gordon inner product 

(<M0 W ') = J dz [<t>Z{oj' + ivd z )4> w i - cj>^, + ivd z )cj)u,] ■ (91) 

The simplified equation (90) possesses a two dimensional conformal invariance. This is 
reflected in the orthonormalized solutions, which are particularly simple: 

<^,± = ^== exp(iw J dzg±{z)) , (92) 

where 

*M=«$±-c-' (93) 
In anticipation to the situation where we shall use non-trivial dispersion relations, let us 
relate g to the Hamilton Jacobi (HJ) action S, a solution of g^d^SduS = 0, see Eq. 
(26). As we shall see classical considerations phrased in the HJ language are useful to 
understand what is modified and what is not when using a non-trivial dispersion relation. 

When considering solutions at fixed energy ui, the HJ action obeys the equivalent of 
Eq. (90): 

(co-vd z S) 2 = c{d z S) 2 . (94) 

Hence the relation between g and S is 

S±{z) =uj J dzg±(z). (95) 

These two actions describe outgoing + and infalling — null radial geodesies. Explicitly 
these are obtained from S by considering W = S—uit and requiring that W be independent 



3 To make contact with Eqs. (36), note that x = ct±cj dzg^. Therefore e '"'^ = 
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of u>. The geodesies are obtained as the locus of stationary phase upon building wave 
packets in ui (the geometrical approximation). 

It is of interest to explicitly consider the behavior of stationary modes in the near 
horizon region, wherein v can be accurately approximated by Eq. (37). In the linear 
regime, 

S + (z) = ^-\n{z-z H ). (96) 

This logarithmic behavior weighted by the frequency in the unit of the surface gravity is 
the key feature characterizing the phase of (outgoing) modes when approaching a (future) 
horizon. On one hand, at the classical level, cL W = gives t — to = \n(z — zh)/ck, i.e., the 
exponential law of Eq. (3) and Eq. (39). On the other hand, this logarithm implies that 
the asymptotic outgoing modes of Eq. (92) are singular on the horizon and only defined 
for z > zh- 

Hence to obtain a complete set of outgoing modes which covers both sides of the 
horizon, one should introduce a new set of outgoing modes living inside the horizon. We 
shall call them trapped modes. A convenient basis of these modes is provided by the 
mirror reflected solutions of the outside solutions: 

<f>Z,tra P = 8(zh - z) x 1 exp(i— \n(z H - z)) , (97) 
**J2lu2-k en 

where 9 is the Heaviside function. These modes have been complex conjugated, because 
they have a negative KG norm when u) is positive. In the Dirac language they would de- 
scribe holes, i.e., in second quantization, they must be associated with creation operators. 
The origin of the sign flip of the KG inner product can be seen from Eq. (91): inside the 
horizon, the sign of p = —id z flips when acting on \n\z — zh\- 

One should also notice that the logarithmic behavior of modes implies that the local 
value of the longitudinal momentum p behaves as 

p = d z S+ = — — !— = Po e- CK(t - to) . (98) 

CK Z — ZH 

The growth of p as 1/ (z — zh ) coincides with the growth of frequencies of stationary out- 
going modes of Eq. (92) as measured in the fluid's atoms frame. Indeed this z dependent 
frequency is 

Q(z) = oj-v(z)p = cp=- - — . (99) 

k z — Zh 

The first equality gives the relationship between the stationary laboratory frequency and 
that in the comoving frame. The second equality follows from the dispersion relation Eq. 
(94). It is interesting to point out that Eq. (99) also appears in the context of gravitational 
BHs, see Eq. (2). In that case it gives the relationship between the stationary asymptotic 
frequency ui and the frequency measured by a free falling observer in the near horizon 
region. With these remarks we have completed the list of classical properties governing 
outgoing propagation in this region. 

Unlike what we just obtained for the outgoing modes, the infalling modes, the — 
solutions of Eq. (92), are regular and cover both sides of the horizon. 
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In summary, having a complete set of stationary modes, which cover both sides, we 
can decompose the field operator as 



■ = 4>- + 4>+ 



f 

Jo 



dcu {a w - 4>u- e ' wt + h.c.) 



)o 

4>+ = 4>outg + (ptrap 
roc 

4>out g = duj (a w , ut g 4>u>,+ e~'"* + h.c.) 
Jo 

<f>tra P = / du) (al trap (j>t ttrap e~ 1 "" + h.c.) . (100) 
Jo 

In the first line we have decomposed the operator into its infalling — and outgoing part +. 
In the third line, we have split the + part into the operators which live outside and inside 
the horizon respectively. The three families of destruction operators a W} -, a w , utg, aw,trap 
annihilate the out vacuum at late time. The Hawking flux will be generated by the action 
of creation operators ai, ou tg- 

4.2.3 The vacuum state resulting from a collapse 

To obtain the Hawking flux, we now turn to the second aspect of the calculation, namely 
the identification of the state of the field which results from the formation of a sonic black 
hole starting in the late past with a subsonic fluid configuration when initially no quanta 
were present. This state is usually called the in vacuum, a denomination we shall use. 

To be able to relate this state to states with a definite particle content built with 
the al operators acting on the out vacuum, we need to consider a metric which describes 
a black hole formation. Since we shall be interested only in the late time behavior of 
Hawking radiation, we can choose the simplest model to describe it. Indeed, any other 
regular way to form the same asymptotic metric will differ from this one by transient 
effects governing the outset of Hawking radiation. The regularity of the process ensures 
that these transients fade out exponentially fast in time in the units of 1/ck. 

The pedagogical model we shall use describes the formation of the sonic hole by a 
sudden transition between subsonic and supersonic motion. For simplicity we assume that 
the transition occurs at a given time to. We assume that initially the fluid is at rest (i.e. 
v = 0), but we could have equally attributed it a finite velocity without affecting the 
result. Then an instantaneous transition causes the formation of the sonic hole in a de 
Laval nozzle. So for t < to v — 0, whereas for t > to, v(z) is given by the stationary sonic 
black hole solution of Eq. (44). A spacetime diagram of the model is shown in Fig. 4. 

The acoustic spacetime is flat for t < to. Hence vacuum positive frequency outgoing 
modes are simply plane waves 

0n = e e ' = e »«' . (I0l) 

V2ty20. V2it2Q 
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Figure 4: Spacetime diagram of an acoustic black hole produced by an in- 
stantaneous shock at t = to. Before the shock the null geodesies propagate 
like those in flat Minkowskian spacetime, but after they are no longer de- 
scribed as straight lines since the spacetime has been modified by the shock 
and is described by the curved acoustic metric (22) (the vertical solid line 
represents the sonic horizon). At late times, the outgoing null geodesies 
coincide with those of Figures 1 and 3. 
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To find Hawking radiation one needs simply to paste these modes to the stationary outgo- 
ing modes <ft u ,+ , 4>Z,tra P of Eqs. (100) which live on either side of the horizon. This is done 
in two steps. First by constructing the in modes after the shock and then by computing 
the overlap of the latter with the stationary modes. 

The result is the following. (The proof is given afterwards.) In the future of the shock, 
unlike the outgoing modes with fixed frequency w, in modes are regular across the horizon 
z = Zh- More precisely, at fixed t they are still of the form, see Eq. (101), 

<fa ~ e lQzh , (102) 

where 7 > is a constant, up to non- universal correction terms which vanish when z — > Zh- 
These describe the outset of HR in the particular formation model one is considering. 

From Eq. (102), we learn that in modes contain only positive momentum p = — id z 
in the near horizon region. This extreme simplicity allows to characterize the in vacuum 
in several instructive ways. Moreover these only refer to the near horizon stationary 
behavior of in modes and thus no longer refer to the region before the shock. Nevertheless 
to consider the latter was necessary to derive the regular properties of in modes. At this 
point it is interesting to mention that this will no longer be the case when considering 
non-linear dispersion relations. In these cases, it is possible to derive HR without any 
reference to the formation of the BH. 

[Proof of Eq. (102). For the first step one needs the relation between laboratory t, z 
coordinates and the null coordinates x~ in both regions. We have, for t < to: 

xTn = ct - z (103) 
whereas for t > to, in the near horizon region we have: 

x~ ut = ct-c [ -^-~ct--\n(z-z h ) + 0(z-z h ). (104) 

J C + V K 

Evaluated along the shock these become 

•^in ^0 z 1 

x~ ut ~ cto-- ln(z - z h ) + 0(z - z h ) , 

HI 

= — — ln(x H — xj n ) + const , (105) 

where X- = cto — Zh is the value of x~ n at which the horizon forms. In the last equality 
we have used Eq. (103) to replace z by x~ n . 

Eq. (105) gives the desired relation x~ ut (x~ n ) which allows to extend the null geodesies 
x^ n = const, in the future of the shock. One has 

x H - Xi„ oc e~ KX °^ oc (z - z h ) e~ CKt , (106) 

up to irrelevant constants. This relation holds in the near horizon limit whatever (regular) 
model for the sonic hole formation is used. At fixed t, this implies that an initial mode 
^ e -^, r n/ c behaves as specified in Eq. (102). QED.] 
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4.2.4 Near horizon in vacuum and HR 



To efficiently characterize the in vacuum and derive the properties of HR, it is particularly 
appropriate to work in the p-representation, even though a priori this might seem rather 
abstract. In fact, these settings are also appropriate to re-derive HR using a non-trivial 
dispersion relation, and to understand the origin of its robustness [19]. 
The Fourier transform of a mode 

<p u (p) = ( dz-^e- lpz <t>^z) (107) 

J V27T 



obeys 

(u — pv)(to — vp)<p w = c 2 p 2 ip u , , (108) 

where v = v(z) is a differential operator acting on its right. In the near horizon region, 
for zh = 0, one has v = — c + icK,d p (see Eq. (37)). Thus Eq. (108) is a second order 
differential equation in d p . In this region, the general solution describing outgoing modes 
is thus 

<p u (p) = A„ e(p)p- i ^ CK - 1 + B w 0{-p) (_p)-W««-i . (109) 
Using Eq. (91), one finds that the Klein-Gordon norm of this solution is 



<<A,|¥V> = J dp2p V l<p u , = (\A LU \ 2 - \B^\ 2 ) (2ty2ck) S(cu - lo') . 



(110) 



From this norm, we learn that positive norm (outgoing) modes are principally made with 
positive momentum. 

The in modes are the simplest ones: they are made exclusively of positive momentum, 
i.e., Bui = for both signs of to, in agreement to what was find in Eq. (102). The 
normalized in modes of eigen-frequency u) are thus 

^, in ( JP ) = e(jp)^2=p- i -'^-\ (in) 

V27r2c/t 

It is now straightforward to get the coefficients of the Bogoliubov transformation 
relating these modes to the singular asymptotic modes. It suffices to inverse Fourier 
transform <fioj,in(p)- Let us proceed in two steps: 



f°° 1 

<?Wn(z) = / dp — 

Jo V2 



' —iuj/c 



V 



V2tt V2tt2ck 
I ^^j^/te^^/". ( 112 ) 
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In the second line, e is a positive infinitesimal number. This equation tells us that the in 
modes are analytic in the upper half z plane, a condition which directly follows from the 
fact that only positive p entered in the Fourier transform. This is crucial: it tells us how 
to connect the singular outside and trapped modes of Eqs. (100) by specifying that the 
cut of lnz is such that the argument of negative z is +in. This fixes the relative weight of 
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the modes living on either side 
1 



0w,in(z) = 



( )e ' \8(z)z ' +e ' 6{—z)(—z) ' 

CK L J 



27tV2ck ck 
^/2-kck/lo ck 

UJ <t>u,,+ + fiu> 4^, trap ■ 



+ C 



0{-z)<t>l,tr 



(113) 



In the second line we have used the normalized outgoing modes of Eqs. (92, 97). The 
Bogoliubov coefficients are, by definition, the coefficients of these modes. They obey 



!/V<M =e" 



(114) 



Together with the unitary relation \aj\ — = 1, one gets Hawking's result: the 

emission rate of outgoing quanta per unit asymptotic time and per unit frequency is given 
by 6 

n w = |/3 w | 2 = (e 2 ^ /cK -l)- 1 , (115) 

that is, a Planck distribution with a temperature given by Th = Hck/2'kkb. 

There is an interesting way to re-write Eq. (113) which is useful to understand the 
robustness of HR. It concerns the decomposition of asymptotic modes in terms of in 
modes. The very fact that the asymptotic stationary modes vanish inside the horizon, for 
z < 0, implies that they are analytic in the lower half p plane. Indeed one has 



</V + (p) 
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dz e 



^r(^ + i )e 

2tvV2^ ck 



(116) 



{(p + ie)-^— 1 } , 



-^—re- + i)c" 

2ttV2lu ck 



,/2c 



x \e{p)p- w/CK - x - e-^ /CK 9{- P )(-p)-^ /CK - ll \ [ . 

From this equation one gets the Bogoliubov coefficients of the inverse transformation. 
They also obey Eq. (114). Hence HR can be equally deduced from the above relation. In 
fact Hawking originally derived black hole radiance in a similar way [2]. He propagated 
backwards in time asymptotic outgoing modes and decomposed them in terms of in modes 
to get the Bogoliubov coefficients. 



6 For the unfamiliar reader, we note that, for free fields, the relationship between in 
vacuum expectation values of observables and Bogoliubov coefficients is universal. Of 
central interest in stationary cases is the mean particle density flux n w , i.e. per unit 
time and per unit frequency, in the in vacuum. It is given by the weight of the delta 
in (0, in|oj,a w /|0, in) = n w 5(u — u>'). In terms of the Bogoliubov coefficients, one has 
n u = \/3 u \ ■ For a detailed derivation of these relations, see [42]. 
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We have now all the tools to address the calculation of HR with the relativistic 
equation replaced by a non- linear Dispersion Relation (DR in the following). Before 
that, let us list the essential properties of modes which fix the properties of HR, namely 
thermality and stationarity. 

A. The modes characterizing in vacuum are regular on the horizon and contain only 
positive momentum p at fixed time, see Eqs. (102, 112). 

B. The stationary modes associated with outgoing asymptotic quanta are singular on 
the horizon and vanish inside it. The singularity of a mode of frequency to is governed 
by a logarithm of z — zh weighted by lj/ck, with k = c 2 k the surface gravity. Since these 
modes vanish inside the horizon, in the p-representation, they are analytic in the lower 
p-plane, see second line of Eq. (116). 

C. This analytical behavior fixes the ratio of Bogoliubov coefficients relating the 
asymptotic modes to the regular modes characterizing the in vacuum, see third line of 
Eq. (116) and Eq. (114). 

These three properties will be essentially unmodified when replacing the relativistic 
equation (88), by that associated with a non-linear DR. This explains why the properties 
of Hawking radiation are robust against introducing a non-linear DR in the UV sector, 
i.e. for frequencies much higher than ck. 

The usefulness of the p-representation can now be appreciated. In the near hori- 
zon region, the mode equation (108) in Fourier transform is second order, as in the z- 
representation. However, the symmetry between the z and p representations is broken 
upon introducing a non-linear DR: in the z-representation the order of the differential 
equation will depend on the choice of the DR, whereas in the p representation, the differ- 
ential equation will stay second order. This allows to perform a general WKB analysis to 
interpret the corrections and to evaluate them. Secondly, for all DRs, the condition for 
being in the in vacuum is still given by imposing that the modes only contain positive 
momentum. 

The important consequence of these two facts, is that the only source of possible modifi- 
cations of HR arising from the near horizon region is backscattering in p-space. 

5 Short distance effects and dispersion 

The first idea is to replace the relativistic relation, which reads 

Q 2 — c 2 p 2 = c p ■ p, (H7) 
in fiat spacetime, by a non-linear dispersion relation 

n 2 = c 2 F 2 (p 2 ), (118) 

which breaks Lorentz invariance, but which preserves the other symmetries of Special 
Relativity, including isotropy. The non-linear behavior of F introduces a UV scale p c 
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which weights the first non-linear term. Generically one has only two cases 



fV)= p 2 ±4 + o(/). 

Pc 



(119) 



The + case gives rise to superluminous propagation (faster than light or sound), whereas 
the — case leads to sub-luminous propagation. This follows from the relation giving the 
group velocity: v g (p) = cdF/dp. 7 

The introduction of the non-linear function F can be conceived from two different 
point of views. First, from the acoustic point of view, when considering the propagation 
of collective excitations in condensed matter, there is dispersion at short wavelength. In 
other words, the Eulerian action of Eq. (5) governs only the evolution of long wavelengths 
configurations. For example, superluminous dispersion is found in dilute Bose gas conden- 
sate [45], whereas the subluminous case is found in 4 He, see [17, 11]. The second point 
of view is more conjectural and related to the short distance behavior of space-time itself. 
Certain authors consider that Lorentz invariance might only be an approximate symmetry 
in nature [17, 46], and will therefore be violated in the UV, say near the Planck scale. 

The second idea is to introduce this modified dispersion relation in a BH spacetime. 
Since the modified DR breaks Lorentz invariance locally, one must choose a particular 
frame. The choice of this frame differs according to the point of view one adopts. When 
one considers sound propagation, there is no ambiguity: the preferred frame is that of the 
atoms of the fluid [17, 18]. In other words f2 in Eq. (118) is the frequency measured in 
the frame of the fluid's atoms. The alternative point related to Quantum Gravity is again 
more conjectural. It is based on the double assumption that the Quantum Gravitational 
effects, which have been so far neglected, will preserve the regularity of the in vacuum 
across the horizon, but will also effectively introduce a non-linear DR. In this second point 
of view, one should choose what is the preferred frame, or better, the preferred class of 
frames, since two frames in a same class will give essentially the same results. The choice 
usually made is that the prefered frame is a freely falling frame. This amounts to introduce 
the modified DR in the PG coordinate systems, see Eq. (199) in Appendix C. 

No matter which point of view is chosen, one deals with the following action. Under 
the same simplifying hypothesis of the previous section, the modified action related to the 
DR of Eq. (118) is 



There is a also third situation, which corresponds to dissipation. This case has not 
been fully analyzed even though it has been represented in [19] and schematically described 
in [44]. The main difference between dispersion and dissipation is that the dissipative case 
is not autonomous: the environment giving rise to the dissipation should be explicitely 
included in order to get a unitary description of this case. We hope to describe this case 
in the near future. 



S = ^ 




(120) 
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Since the modification concerns only the spatial part of the action, the inner product and 
the normalization condition, Eq. (91), are both left unchanged. 
The modified mode equation, replacing Eq. (90), becomes 

(w + id z v)(u + ivd z )<t> w = cF 2 {-d 2 z )4>^ ■ (121) 

To study the near horizon behavior of mode propagation, we return to the p-representation. 
The modified equation in Fourier transform reads 

(u> — pv)(tu — vp)(p W} F = c 2 F 2 (p 2 ) <p W} F (122) 

where v — — c + icK,d p , as in Eq. (108). 

The structure of this equation guarantees that for all F the solutions can be written 
as a product 

^,F = ^xe- lp/K X F, (123) 

wherein the first factor is the solution of the relativistic equation, hence F independent, 
whereas the second factor is uj independent. The function xf obeys the simple second 
order differential equation: 

«V d 2 X f = F 2 X f. (124) 

Moreover, when <fiu>,F is conventionally normalized with the KG product (which is F 
independent), see Eq. (110), xf must have the following unit Wronskian 

-^{XF-idpXF ~ XFidpXF) = ±1- (125) 

The + sign applies to outgoing modes, whereas the — sign applies to infalling modes. 

The equations (124) and (125) fully govern the modifications of HR which come from 
the near horizon region. Of course there will be other corrections which arise from the 
outer region wherein the gradient of the velocity goes down to zero, see e.g. [47]. However 
these modifications will affect the propagation of modes irrespectively of HR, as the 'grey 
body factor' does it in the Schwarzschild metric in the 3-dimensional case. In what follows 
we shall only discuss the (intrinsic) modifications of HR which arise from the near horizon 
region. 

To present the nature of these modifications, we shall split the presentation into three 
parts. First we discuss the mathematical properties of xf and its WKB approximation. 
Then we present the calculation and the interpretation of the corrections to this approxi- 
mation. At the end only we argue why the asymptotic conditions to determine the modes 
which describe in vacuum and asymptotic quanta are both unaffected by the non-linear 
behavior of F. 

5.1 Mathematical properties and WKB approximation 

Let us consider the exact solution of Eq. (124) which for p —* oo is purely 'forward', 
i.e. it is a running wave which possesses a unit positive Wronskian of Eq. (125). When 
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evaluated for moderate values of p, i.e. p/p c <S 1, the exact solution will be of the form 

— A ~WKB , R WKB,* (T>a\ 

where Xf K B is the following unit Wronskian WKB solution 

The conservation of the Wronskian guarantees that 

\A F \ 2 - \B F \ 2 = 1 . (128) 

In the WKB approximation, Bf = 0: no backscattered wave (in p) is engendered in 
the journey from p = oo down to p <C p c . In this regime, the norm of Af stays constant: 
|yljr(p)| = |^4f(p = oo)| = 1. Hence, there will be no modification of HR in this WKB 
regime. Indeed the ratio of the Bogoliubov coefficients, determining Hawking temperature, 
is extracted from the relative weight of the positive and negative norm modes and not from 
the local behavior of the modes. This local behavior is in fact affected by the non-linear 
behavior of F: up to a constant phase one has 

e -^ K X F KB = l + 0(p/p c ) 2 . (129) 

The important point is that the correction terms in (p/p c ) 2 should not affect the Bogoli- 
ubov coefficients of Eq. (113). 

5.2 Backscattering and modifications of HR 

The modifications to HR come from the global corrections to the WKB approximation, 
i.e. the amplitude Bf. Indeed, there are two kinds of modifications to Xf KB '■ first there 
are local (in p) modifications which follow from the fact that Xf /KB does not exactly obey 
Eq. (124) at that p. These corrections should not affect HR, see the above discussion. 
Secondly there are global (accumulated) errors which result from the fact that Xf KB did 
not exactly obey Eq. (124) from p = oo down to values p ~ k. This second type is 
responsable for Bp 7^ in Eq. (126). 

To evaluate Bf is not an easy task because it is an accumulated error, see e.g. [48]. 
What can be said is that the local source of Bf is 

This dimensionless quantity is smaller than n/p c . From this we learn that scale separation, 
i.e. n/pc <JC 1, is, in full generality, a sufficient condition for having small corrections to 
HR. 

Moreover, since the final value of Bf (i.e. for p of the order of k) is expressed as 
an oscillatory integral of cr(p), the integral will generically be smaller than the integrand 
(= a). To illustrate this we consider an exactly solvable case: the quartic dispersion 
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F 2 = p 2 +p 4 /p 2 . In this case, Eq. (124) describes an up-side down harmonic oscillator. 
Then using the equations of Sections 19.17 to 19.22 of [49] one gets 

B F = e -^ /R , (131) 

which is indeed much smaller than a ~ n/p c . 

The lesson is that the (intrinsic) corrections to HR originate from the backscattered 
amplitude and not from the local corrections to Xf KB ■ The backscattered amplitude is 
small whenever the UV scale governing dispersion is much larger than the surface gravity. 

A more subtle consequence of Eq. (126) is that Af, the amplitude of the forward 
wave, has increased when Bf =fc 0. Hence the corrections to HR are related to another 
pair creation process. The new Bogoliubov transformation mixes outgoing and infalling 
modes, unlike what was found in Eqs. (113). To our knowledge this has not been described 
in the literature. We hope to report on this in the near future. 



5.3 The identification of in and out modes 

It remains to verify that the non-linear behavior of F does not affect the relationship 
between the asymptotic physical states and the mathematical properties of the modes, 
solutions of Eq. (121). 

Let us start with in vacuum. We claim that the initial characterization of the in 
vacuum is still that the in modes only contain positive p. Let us explain why in the light 
of the properties of the classical propagation in the near horizon region. 

The factorization of the modes in the p-representation is recovered in the HJ formalism 
in the following way. The modified Legendre transformed HJ action, Wf(p) = Sf(x) —px, 
obeys the equivalent of Eq. (121), 

Lj+p{c + CKd p W F )=cF{p 2 ). (132) 

Its solution is a sum of the relativistic action plus a modified term independent of u> 

Wf{p) = ^ 1qp+ r dp > F ^)-p\ (133) 

CK J Kp 

This tells us that, irrespectively of the choice of F, the trajectory in p-space is, in the near 
horizon, the same as that of the relativistic case, namely Eq. (98). So, even though the 
trajectories in spacetime do depend on F, in p space they always obey the exponential k 
decay. 

It is important to analyze the propagation in z space to see that the exponential 
blue-shift (backwards in time) is cut off by p c , the UV scale brought in by F. Eq. (98) 
translates in z space in the following manner: the difference in z between the classical 
trajectories with positive and negative energy u> always obeys 

zfA*) - «F,-u,(t) = AzfA*) = AzF^(t) = — e CKt . (134) 
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Therefore the characteristics of the modified dispersion relation are still turned apart by 
the surface gravity, as they were in the relativistic case. However, the important difference 
is that this stretching lasts only a finite time for non-linear dispersion, either super- or 
sub-luminous. To prove this we consider the motion of the "center of mass" of a couple 
of trajectories uj, —uj. It is given by 



It is independent of w, as was the behavior of Xf- In the relativistic case, F = p for all p, 
one recovers z = 0, the locus of the static event horizon. Prom Eq. (135) we learn that 
the couple leaves the near horizon region, i.e. \zf\ = for p ~ p c . Therefore Eq. (98) 
applies only until p ~ p c . This is the main effect of the non-linear behavior of F: when the 
critical momentum is reached, the center of mass starts to drift and this brings the pair 
of characteristics outside the near horizon region, where the gradient of v progressively 
vanishes. 

Imagine indeed that the gradient of the flow smoothly decreases as one leaves the near 
horizon region. Hence the acoustic metric will become flat and p will tend to a constant. 
In this steady regime, far away from the horizon, positive p characterizes in vacuum for 
outgoing modes, irrespectively of the dispersion relation F. Hence this proves that the 
in vacuum is still asymptotically characterized by positive p for p ~ p c 3> k. It should 
be emphasized that we did not refer to the formation of the BH to characterize the in 
vacuum. This is no longer necessary because the blueshifting effect stops after a finite 
time and the field configurations at the origin of HR leave the horizon region and reach a 
flat region where a usual characterization of vacuum by plane waves applies. Instead, in 
absence of dispersion, the blueshifting effect only stops when the horizon forms. 

The characterization of outgoing asymptotic quanta is clearly unaffected by F because 
it is done in a regime where p ~ k <§C p c . The recovery of the unperturbed relativistic 
behavior is illustrated by the fact that when p <C p c , zf oc p 2 /pi — > 0, i.e. zf tends to the 
position of the horizon, as expected. 

From this analysis, we have learned that the characterization of both in modes and 
asymptotic out modes is unaffected by the introduction of F. Hence the only source of 
modification of HR comes from Bf 7^ in Eq. (126). 

5.4 The lessons 

Since several lessons have been drawn by the analysis of HR in the presence of a non-linear 
DR it is of value to present them synthetically. 

I. HR is robust. 

The properties of HR are robust in the sense that stationarity is fully preserved, whereas 
the modifications of the spectrum are smaller or equal to n/p c <C 1. 

II. The origin of the robustness of HR. 

It is also of value to specify why HR is so robust. In the near horizon region, the modified 



zf (t) + z F , 
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= 2F,a>=0 = 



F( P )-P 



(135) 
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modes can be expressed as the relativistic modes times a running wave function xf which is 
independent of the frequency w, see Eq. (123). Hence when there is no backscattering and 
the three properties (A, B and C presented at the end of Section 4) of the relativistic 
modes are unmodified, no modifications of the Bogoliubov coefficients appear. What 
governs the backscattered amplitude is a non-adiabatic effect, like the Landau-Zener effect, 
which can even be exponentially small, see Eq.(131). 

III. The exceptional feature of the relativistic propagation. 

From Eq. (134) we learn that the endless stretching found in the relativistic case is 
an exceptional feature resulting from the absence of any non-linearity in the DR. This 
maintains the pair of characteristics centered for all times on the event horizon z = 
0. Instead, for all non-linear DR, the departure from to = cp as p increases induces a 
drift which sends the pair outside the near horizon region where the gradient of velocity 
vanishes, hence where the blueshifting effect ceases. 

IV. Hints for Quantum Gravity. 

This last remark is an invitation to inquire about the physical validity of the relativis- 
tic equation Q. 2 = c 2 p 2 to characterize field propagation near a horizon. This equation 
results from the hypothesis that the bare actions are locally Lorentz invariant, but also 
from a linearized treatment. Without questioning the validity of Lorentz invariance at 
a fundamental level, it is possible that radiative corrections evaluated in the non-trivial 
near horizon metric will modify 1-loop Green functions in such a way that the effective 
law governing their behavior will be characterized by a new UV scale which shall act like 
Pa did it in the non-linear DR we used. A first attempt to show that this might be the case 
has been presented in [24]. We refer to [51, 52] for a more fundamental but less ambitious 
approach to derive, from radiative corrections, effective dispersion relations breaking the 
(local) Lorentz invariance possessed by linearized field propagation. 

6 The phonons backreaction 

In the mean, the backreaction of the linearized phonons on the underlying fluid can be 
discussed using the sound equations that govern the classical backreaction [20, 21], i.e. 
Eqs. (61, 62) (see also [50]) 

ps+V(p s « s )+VV ,T 5 5 <Tf)=0, (136) 
^B + \vi+ l x{p B )-\^f{T)=Q, (137) 

where now (T p „) are the quantum expectation values of the pseudo energy momentum 
tensor operator of the phonon field. Because of the analogy, the (T M „) coincide with the 
expectation values of the stress tensor for a minimally coupled scalar field propagating in 
an effective acoustic spacetime whose metric is g^(pB,VB) of Eq. (56). 
In the case of our interest, these expectation values are to be taken in the in vacuum, 
in which in the remote past prior to the time dependent formation of the sonic hole, the 
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phonon field is in its vacuum state. 

The quantum version of the Euler equation, namely 

'v B + {v B ■ V)v B + V/x(p B ) - {^f{T)) = (138) 

is the fluid analogue of the semiclassical Einstein equations (see Appendix D) in gravity. 

The analogue of the continuity equation (61) is the covariant conservation of the stress 
tensor in Einstein's General Relativity. 

It is instructive to compare the role played by the stress tensor in determining the dynam- 
ical evolution of the system. Unlike semiclassical Einstein gravity, where all components 
of {Ty, v ) enter the dynamical equations (see Appendix D, Eq. (210)), in hydrodynamics 
the dynamical backreaction is caused just by the gradient of the trace density. 
As they stand, the expectation values of the minimally coupled scalar field ipi entering 
the backreaction equations (61, 62) are formally ultraviolet divergent. 8 
Here we are touching a very delicate point. In performing the integrals over the frequen- 
cies involved in the expectation values of one inevitably enters a regime where 
the wavelength of the modes becomes comparable to the intermolecular separation. The 
approximation of the fluid as a continuous medium breaks down. High frequency modes 
see the elementary constituents of the fluid and the dispersion relation for these modes 
deviates significantly from the relativistic one (iv = c\p\) associated to the wave equation 
(20) as we have discussed in section 5. We make the assumption that the contribution of 
these modes to the expectation values is negligible. We shall give later on arguments in 
favor of this. Our working hypothesis is therefore that the wave equation (20) holds for 
all frequencies. Furthermore, we assume that the divergent expectation values are made 
finite by adopting the standard covariant renormalization scheme in curved space (see for 
example [36]), which, as such, preserves at the quantum level the covariant conservation 
law of the classical pseudo energy momentum tensor (i.e. {T^) :lJ = 0). Within this frame- 
work one can eliminate the divergent part of (T^J). 

Unfortunately the finite part of (T p „) for an arbitrary metric is not known and therefore 
the actual acoustic metric, solution of the backreaction equations, can not be determined. 

6.1 The dimensional reduction model 

Lower dimensional models are very useful in physics as they allow explicit analytic so- 
lutions of the dynamical equations to be found, providing at least a qualitative insight 
of the behavior of real four dimensional systems. This attitude has been largely used in 
studying quantum effects in black holes [54, 56]. The basic assumption consists in re- 
placing the four dimensional quantum stress tensor by a two dimensional one associated 
to a two dimensional minimally coupled scalar field propagating in the two dimensional 
section of the physical four dimensional metric. As we know, this approximation neglects 

8 Infrared divergences can be simply eliminated by introducing a mass cutoff [53]. 
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the backscattering of the quantum field caused by the potential barrier. 
We follow the same attitude for a qualitative insight in the quantum backreaction occur- 
ring in a sonic black hole [21, 20]. 

We shall focus on a supersonic flow through a de Laval nozzle, which can be considered 
as one dimensional as discussed in subsection 2.6. Under these assumptions, we can write 
the classical equations of motion for the fluid as: 

Ap + d z (Apv) = (139) 
a{^+\v 2 +^p)^ =0 (140) 

where v = d z ip and we have also assumed the time independence of the area A of the 
transverse section of the nozzle. Proceeding as in section 2.7, we separate the perturbations 
from the mean flow. Linear perturbations, i.e. sound waves or phonons at the quantum 
level, satisfy the wave equation 



d a [^^/^9 Bab d b ^yO (141) 

where 



g °V> = -PZl *—* "° | (142) 

c y v B -l 

is the two dimensional (t, z) section of the metric (56) and a,b — t, z. Note the 
invariance of Eq. (141) under Weyl transformations. 

The two dimensional backreaction equations for the mean flow Pb,vb read: 



Ap B + d z (Ap B v B ) + d z 



, (143) 



A Ub + \v% + Mp s )J - i(T (2) ) = 0. (144) 
Here T™ = ^(d a ^id b ^ - \g (2)cd d c ^d d ^g%). 

As seen from Eq. (141) the quantum field ipi behaves as a massless scalar field which is 
coupled not only to the two dimensional metric but also to the "dilaton" This 

causes backscattering of the modes. The Polyakov approximation consists in neglecting 
the backscattering. Eq. (141) is then approximated as 

□ (2) t/>i = , (145) 

i.e. ipi is treated as a two dimensional minimally coupled massless scalar field. The reason 
of this simplification is that even for the "exact" two dimensional theory of Eq. (141) the 
exact (Tab) is not known. Only its trace is known being completely anomalous [55]. Given 
this assumption, the quantum expectation values entering Eqs. (143, 144) can therefore 
be evaluated using the methods developed in Appendix E. 

We shall use Eqs. (143, 144) to find the first order in h corrections to the classical fluid 
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flow in a de Laval nozzle which corresponds to a sonic black hole as described in subsection 
2.5. To this end one writes the dynamical variables pB and vb as 



Pb = Pd + p q , (146) 
vb = v c i + v q (147) 

where p c i, v c i satisfy the classical equations (139, 140) and p q , v q are the corrections of 
order h induced by the quantum fluctuations of the phonon field encoded in (T^ 2 }). 
We can therefore write the backreaction equations (143, 144) as 



Ap q + dz [A(p q Vd + PclVq)] - 8z 



-I((^))+, d (Tif)» 



, (148) 



c 



2 



< T (2)> 



A(lP q + V ci Vq + p q ) - ±— -!- = . (149) 

pd 2, 

Up to the order O(h), the {T^ 2 }} are computed on the classical solution p c i, v c i, which are 
given in Eqs. (44, 47), namely 

A ^JLA H e^- c2 ^ 2c2 , (150) 

Pd = poe"^ /2c2 . (151) 

This configuration, as discussed in subsection 2.6, should be regarded as the final config- 
uration of the fluid resulting from the time dependent formation of the sonic hole. 
We will be mainly interested in the region near the sonic horizon z = where the profile 
of the nozzle can be approximated as 

A = A H +Pz 2 . (152) 

For the sake of simplicity, the suffix "cl" will be omitted in the following. 
Standard covariant renormalization in curved space allows, under the previous assump- 
tion, the expectation values of the phonon fields in Eqs. (143, 144) to be evaluated (see 
Appendix E). For the trace term we have simply (trace anomaly) [57] 

<T (2) ) = ^R (2) , (153) 
whereas in the continuity equation (143) we can write 

where we have used the null (retarded and advanced) coordinates (36) 

X±=C { t± l^)- (155) 
In our double null coordinates the components of the stress tensor are (see for example 
[36]) 

(T±±) = —C^CT^l 2 + t± , (156) 
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The conformal factor C of the two dimensional metric (142) reads 

C=£(c 2 -, 2 ) (157) 

and t± are functions of x ± respectively, which depend on the choice of the quantum state 
in which the expectation values are to be taken. 

Near the horizon (x~ — > +00) the t± can be approximated with the ones characterizing 
the in state at late time [58], namely no incoming radiation in the past: 

t+ = , (158) 

and thermal radiation of phonons at the Unruh temperature in the future: 

This latter represents the Hawking radiation flowing at late time at z — > 00. 
For our configuration the surface gravity k reads 

k = ck = <?s\4- ■ (160) 
V Ah 

Expanding the classical background near the horizon z = to evaluate the (T^), one 
eventually arrives at the following solution of Eqs. (143, 144) for the quantum corrected 
velocity 

c 2 2 

vb — v + v q = — c + ckz k z + e(bi + c\kz)kI , (161) 

6 

where e = a2 p g_i/2 e ' s the dimensionless expansion parameter and bi = 97/2, ci = 
-3047/15, 7 = kc 2 A h /24ty, and 

Pb = p+ epq = poe~ 1/2 ^1 + kz — k 2 2 2 /6 + e ^-( a + Sitz)tj , (162) 

where a = — 02+CK61, S — —bi — ai+CKC\, with 0,2 — 11897ck/120, &2 = — 1618177ck/120. 
The solution is valid for kz < 1 (near horizon approximation) and cut <C l. So equations 
(161, 162) give just a hint on how the backreaction starts near the sonic horizon. 
The boundary condition chosen to derive Eqs. (161, 162) is that, at some given time (say 
t = 0), the evaporation is switched on starting from the classical configuration (p,v): i.e. 
v q (t = 0,z) = 0, p q (t = 0,z) = 0. 

Inspection of Eq. (161) reveals that the net effect of backreaction is to slow down the fluid 
(v q > for nz <g; 1). This goes with a decrease of the density (p q < 0) in the same limit. 
The same equation allows to determine the evolution of the acoustic horizon, defined as 
the locus where vb = —c. We find 

ZH = -?f, (163) 

i.e. the horizon is moving to the left: the hypersonic region shrinks in size. The coefficient 
61 , which determines the quantum correction to the velocity field, and hence the evolution 
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of the horizon, is proportional to the gradient of the trace calculated at t = 0, z = 0, as 
it can be seen by taking the gradient of Eq. (144) in this limit. This should be compared 
to the gravitational black hole case where the evolution of the horizon is determined by 
the energy flux (M oc (T r t ) in spherical symmetry). While in the latter case Hawking 
radiation occurs at the expense of the gravitational energy of the black hole (the mass 
M), the phonons emission takes away kinetic energy from the fluid. 

Considering the evolution of the acoustic black hole as a sequence of quasi stationary 
configurations one can find the correction to the emission temperature 



= h_dv 
2n dz 



he 



563 3 

-en cAut 



720tt 



(164) 



As time evolves, the temperature of the emitted phonons decreases. 

As discussed in Appendix D, the opposite happens for a gravitational black hole of mass M. 
Being T ~ 1/Af, as the mass decreases because of Hawking evaporation, the temperature 
grows causing the black hole to evaporate completely in a time of order Mq (Mo is the 
initial mass of the hole). Instead, if the black hole possesses a conserved U(l) charge Q, 
for M ~ \Q\ a temperature decrease is associated to the mass evaporation, which stops 
when it is reduced to M = \Q\ (the so called extremal black hole) for which the Hawking 
temperature is zero. This ground state is approached in an infinite time. Acoustic black 
holes seem therefore to behave like near extremal charged black holes. 



6.2 Discussion 

An important question which remains open is to what extent the above results depend on 
the two assumptions we made, namely, first, the use of a relativistic dispersion relation 
which ignores the short distance corrections due to the molecular structure of the fluid, 
and second, the covariant regularization scheme used to subtract ultraviolet divergences. 
These two aspects are deeply connected. For the Polyakov theory that we have used 
in the 2D backreaction equations (143, 144), Jacobson [22] has argued that, within a 
covariant regularization scheme, no significant deviation from the usual expression for 
the trace anomaly and the components of (T^/) are expected if one introduces a cutoff 
at high frequencies. For the hydrodynamical system we have considered, covariance is 
a symmetry of the phonons low energy effective theory only, which is broken at short 
distance. Hence non covariant terms depending on the microscopic physics are expected 
to show up in the effective action, and these are crucial for a correct description of the 
unperturbed quantum vacuum of the fluid. However the expectation values (T^ v ) entering 
the backreaction equations (136, 137) do not represent the energy momentum of the 
fluid quantum vacuum. They describe instead the perturbation of the stationary vacuum 
(whose energy is strictly zero [11]) induced by inhomogeneities and by the time dependent 
formation of the sonic hole which triggers the phonons emission. Here we have assumed 
that these deviations can be computed within the low energy theory. This situation is not 
unusual. Casimir effects are well known examples of vacuum disturbances caused by the 
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presence of boundaries. It happens that the Casimir energy is often (but not always, see 
G. Volovik in Ref. [30]) independent on the microscopic physics and can be calculated 
within the framework of the low energy theory. As shown by Bachmann and Kempf [59], 
deviations caused by the modified dispersion relation are suppressed by positive powers 
of d c /L, where d c is the interatomic distance of the material constituting the walls and 
L the separation between the walls. This happens because, while low frequency modes 
are reflected by the boundaries, for the high energy ones the wall is transparent. They 
produce a divergent contribution to the vacuum energy which is canceled by a proper 
regularization scheme and does not affect the finite result. Similarly, in [60] it is observed 
that, although in a Bose- Einstein condensate quantum fluctuations have a non-linear 
dispersion relation, the Casimir energy is given by the long wavelength modes, which 
therefore behave as massless scalar fields. We have assumed that a similar decoupling 
happens for the acoustic black hole. The check of our hypothesis would require an analysis 
of the quantum system within the microscopic theory which takes into account the time 
dependent non homogeneous formation of the sonic hole. This is for the moment beyond 
computational capability. Anyway, as it has been shown in section 5, modifications of 
the dispersion relations, to take into account short-distance behavior of the high-energy 
modes, basically do not affect the spectrum of the emitted phonons. This is not a proof 
that observables like (T^) are also unaffected by short-distance physics. However, keeping 
in mind what (T MI/ ) does really represent (i.e. deviations in energy and momentum caused 
by the inhomogeneity of the flow), and the indications coming from the Casimir effect, 
one expects the contribution coming from the relativistic low enegy theory to be the most 
significant in driving the backreaction. 



7 Conclusions and open questions 

In section 2 we show how a relativistic wave equation emerges when studying the linear 
perturbations of an Eulerian fluid which possesses a Galilean invariance. As mentioned 
in the footnote 3, one gets a similar wave equation when considering relativistic fluids. 
In both cases, the wave equation defines an "acoustic metric", i.e. a Lorentzian four- 
dimensional metric which is curved when the fluid's flow is non- homogeneous and/or non- 
stationary. Similar wave equations appear in many non-linear systems, such as non-linear 
electrodynamics [61], and di-electric media [62]. It would be interesting to identify what 
are the physical ingredients (such as local equilibrium) necessary to obtain a Lorentzian 
wave equation for linear perturbations. 

We conclude section 2 by analyzing the first backreaction (non-linear) effects, namely 
how the mean fluid's flow is affected by the energy carried by an ensemble of sound 
waves. We show that these effects are governed by the pseudo energy-momentum tensor 
associated with the four-dimensional metric, see Eqs. (61, 62). It would be also interesting 
to see if this result extends to other physical situations wherein linear perturbations obey 
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a Lorentzian wave equation. 

In sections 3 and 4, we review the essential features of quantum field theory in curved 
spacetime, focussing on the case of interest: the stationary curved metric describing a 
(static) black hole geometry. We emphasize that the stationary modes associated with 
asymptotic quanta are singular on the horizon, while the quantum vacuum state is regular 
there. This potential conflict is resolved by the existence of Hawking Radiation, which is 
encoded in the non-trivial Bogoliubov coefficients relating the singular asymptotic modes 
to the regular modes describing the vacuum state, see Eq. (113). We have phrased 
the derivation in the p-representation to prepare the analysis of HR in the presence of 
dispersion, which appears at short distances in materials. 

In section 5 we explain why the properties of HR, namely stationarity and thermality, 
are not affected when introducing a non-linear dispersion relation, in spite of the facts 
that Lorentz invariance is broken at short distances, and that the field propagation is 
dramatically modified in the near horizon region. From this analysis we learn that the 
properties of HR are robust, but also that the near horizon behaviour of the modes in the 
relativistic case is exceptional (unstable) in the following sense. Any non-linearity in the 
dispersion relation will indeed erase their singular behavior because the blue-shifting effect 
obtained in a backwards in time extrapolation stops when the UV scale is reached, see 
Eq. (135). Possible implications of this result for Quantum Gravity are briefly discussed. 

In section 6 we compute the secular effect on the fluid's flow induced bybHawking 
radiation. The expectation values (T^) of the pseudo energy-momentum tensor of the 
phonons drive this backreaction. Their computation depends in principle on the UV be- 
havior of the dispersion relation. However, there are strong indications that the most 
significative contribution to (T p „) comes from the low energy theory (i.e. relativistic dis- 
persion relation). Given this the evaporation of a sonic black hole seems characterized 
by a slowing down of the fluid associated to a decrease of the temperature of the emitted 
phonons. 

In this analysis the starting point was the action (5) and the fields p and ip were considered 
as fundamental. In many condensed matter systems, however, the hydrodynamical action 
represents just an approximation of the real action, e.g. in Bose-Einstein condensates. 
For these systems, the backreaction problem should be addressed at the level of the real 
action. The correction to this more rigorous approach with respect to the hydrodynamical 
one developed so far is matter for future investigations [63]. 

Acknowledgments: The authors thank S. Liberati, M. Rinaldi and M. Visser for useful 
comments and suggestions. S. F. also thanks Centro Enrico Fermi for supporting her 
research. 
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Appendix 



A General Relativity 



General Relativity is a geometric theory of gravity. According to it the gravitational 
interaction manifests itself as curvature of the spacetime manifold. The gravitational field 
is described by ten potentials which are the independent components of the geometric 
tensor g^v, giving the "distance" in spacetime between neighbouring events 

ds = g^vdx^dx" . (165) 

According to the principle of equivalence, the effects of gravity can be locally eliminated 
going to the local inertial coordinate system in which the metric tensor assumes the 
Minkowski form ?? M „ = diag(— 1, 1, 1, 1) and the connection coefficients 

v % = \d al) + £M,m ^ 9v,f}) (166) 

vanish. 

Gravity, however, shows up as relative acceleration between neighbouring points causing 
tidal distortions measured by the Riemann tensor of the spacetime 

R M _ pM _ pM , pM p7 _ pM p7 fl 671 

In absence of gravitation spacetime is flat and the Riemann tensor identically vanishes 
(and vice versa). 

The sources of the gravitational field are the distributions of mass-energy of the matter 
and gauge fields encoded in their energy momentum tensor T^. This determines the 
gravitational potentials, up to coordinate transformations, via the Einstein equations (here 
and in the following we take G = c = 1) 

~ = 8ttT mi , . (168) 

Here 7? M „ is the Ricci tensor 

R^v = R^ a v (169) 

and R the Ricci scalar 

R = g^R„ v . (170) 
Because of the Bianchi identities the stress tensor is covariantly conserved 

t;, v = o , (171) 

where ; means covariant derivative with respect to the metric <? M „. This describes how 
energy and momentum are exchanged between matter-gauge fields and the gravitational 
field. 

In the vacuum, T M „ = and Einstein equations simplify to 

R^ = 0. (172) 
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Once the gravitational field is determined, by Eqs. (168) or outside the matter energy 
distribution by Eqs. (172), one can analyze the motion of test particles. These move along 
the geodesies of the curved spacetime, determined by 

where t is an affine parameter along the geodesic, which, for massive particles, can be 
chosen as proper time. In term of the tangent vector K a — the geodesic equation can 
also be given as 

K a0 K p = . (174) 

Geodesies of interest are either timelike (with tangent vector normalized to K^K^ = — 1) 
in the case of massive particles, or null (K^K^ = 0) in the case of massless particles (for 
example photons). 



B Gravitational waves 



Gravitational waves are curvature ripples rolling across the spacetime. We shall briefly 
outline how to derive the propagation equations and the characteristic features of gravi- 
tational waves. A more complete and detailed discussion can be found in [3]. 
We shall consider gravitational waves propagating through a vacuum spacetime. Let R be 
the typical radius of curvature of the background, A/2-k and A the typical reduced wave- 
length and amplitude of the waves respectively. We assume A« 1 and A/2-k <C R. The 
spacetime metric g^ v is expanded in a "background" g]S) , a smooth part which included 

(175) 



a (B) 4- h 

Hill/ \ H'liv 



the backreaction of the waves, plus perturbation fo M „, i.e 
The Ricci tensor can be expanded as 



( A 1 



where R.\E) is the Ricci tensor constructed from gffl , while 
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Here and in the following indices are raised and lowered with gJZ 1 , h = g a/3 ( B ) h a g and 



(S) 



(176) 



(177) 



(178) 



(s) 



" ; " indicates covariant derivative with respect to g^ 
The vacuum field equations require i? M „ = 0. The linear part in the amplitude A, i.e. 
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Fuj}} , must vanish since the distorting effect of the wave on the background is a non linear 
phenomenon, so R^J cannot be linear in A. Hence 

fl#(V)=0- (179) 

These are the propagation equations for the gravitational waves. 
The rest of 7? M „ can be written as 

B$ + (R$(h)) + O (£\ = (180) 



averaging over several wavelengths. We can rewrite this in a suggestive form 

1 

2 

where 



Rff-^RWgft^&nTW (181) 



rp(gw) _ 

^ " 8 



l - ((R$hh)) -\g$\R*\h))) (182) 

is the effective stress tensor for the gravitational waves. 

The propagation equations (179) can be further simplified in terms of 

hpv = h^ v - ^hgffl . (183) 

Reparametrization invariance can be used to enforce the "Lorenz gauge condition" 

= . (184) 

In this gauge, the vacuum propagation equations for the gravitational waves becomes 

□ V + 2fl£^ft a " = . (185) 

In the geometric optics approximation, one assumes hpv to have an amplitude that 
varies slowly (on scale I <ti R) and a phase which varies rapidly, namely 

Ku = »(A^„ + ..)e ie . (186) 

One can show that the rays (the curves perpendicular to the surface of constant phases) 
are null geodesies, i.e. 

K a K a = , (187) 
K aiP K fi = 0, (188) 

where K a = 6, a . Being K a a null vector, gravitational waves propagate at the velocity of 
light. 

C Black holes 

Since gravity couples to energy-momentum it affects not only the motion of massive par- 
ticles, but also of light. In the presence of a gravitational field, light does not follow a 
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rectilinear motion, but it moves along the null geodesies (173) of the curved spacetime. 
The deflection caused by the gravitational field on the motion of light is very small in 
"weak" gravitational field like the one of our Solar System. However, there are astro- 
physical systems where this deflection is so strong that light cannot escape from some 
compact region of space that becomes causally disconnected from the rest of the Universe. 
This region is called black hole, the most spectacular consequence of Einstein's General 
Relativity. In a black hole, the gravitational field is so strong that light is dragged back 
instead of escaping from it. 9 The simplest example of a black hole is the one based on the 
Schwarzschild solution to the Einstein equations. This is the unique solution of vacuum 
Einstein Eqs. (172) under the hypothesis of spherical symmetry (according to the Birkhoff 
theorem [64]). 

In the so-called Schwarzschild coordinates the line element reads 

ds 2 = -(1 - 2M/r)dt 2 + (1 - 2M/r)~ 1 dr 2 + r 2 (d6 2 + sin 2 6d<p 2 ) . (189) 

This metric describes the gravitational field of a spherical body of mass M. It is asymp- 
totically flat (i.e. <? M „ — > 77,ii/ for r — > co) and has the correct Newtonian limit. The metric 
presents singularities at r = 2M and r = 0. The latter is a curvature singularity (i.e. 
the invariants constructed from the Riemann tensor diverge), which signals a breakdown 
of the theory there. On the other hand, the singularity at r = 2M is just an artifact of 
the coordinates system used and can be eliminated by a coordinate transformation. For 
example, using the advanced or ingoing Eddington-Finkelstein coordinate [65] 

v = t + r* , (190) 

where r* = J (j^W/V) ~ r + 2Af In \ r/2M — 1|, the Schwarzschild line element can be put 
in the following form 

ds 2 = -(1 - 2M/r)dv 2 + 2dvdr + r 2 dQ, 2 , (191) 

where dQ, 2 — d6 2 + sin 2 6dip 2 is the line element of a unit two-dimensional sphere. 
Due to the cross term drdv in these coordinates the surface r = 2M is regular. It is a 
null surface with special properties, as we shall see in short. The radial null geodesies 
(8 = const., ip = const.) of the metric (191) are found by setting ds 2 — 0; one finds 

v = const. , (192) 

which describe the motion of ingoing spherical light fronts, and 

dr 1 

|J = -(l-2M/r) (193) 

for outgoing light. A plot of both families of geodesies is given in Fig. 5. Ingoing ones are 
depicted as straight lines. 
Inspection of Eq. (193) clearly indicates that outgoing rays manage to move outwards 

9 Since nothing can propagate faster than light, the same thing happens to any matter 
field. 
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Event 



Singularity 




Figure 5: Diagram representing the radial geodesies tangent to the light cones 
in the spacetime of a Schwarzschild black hole. On the vertical axis we have 
increasing times, while on the orizontal one we have the radial coordinate 
r (two angular coordinates have been neglected). In r = there is the 
curvature singularity, while the horizon is located at r = 2M. The region 
r < 2M is the gravitational black hole. The light cones are increasingly 
distorted when approaching the horizon. Once the horizon is crossed, they 
all point towards the singularity. This shows why nothing can ever escape 
from the black hole (see also Fig. 6). 
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only for r > 2M (for which $L > 0), whereas, for r < 2M, they move inwards < 0)! 
For r = 2M they are stuck there forever. 

Therefore, light and any matter field cannot escape from the region r < 2M, which 
represents the black hole. The boundary of this region (i.e. r = 2M) is the horizon, 
which can be crossed only from outside to inside. The spacetime given by the line element 
in (191), and whose structure is depicted in Fig. 5, describes, according to relativistic 
astrophysics, the gravitational field exterior to a sufficiently massive spherical star (M ~ 
3M sun ) at the end of its thermonuclear evolution, when the internal pressure is no longer 
able to prevent the star collapsing because of its own gravity [66]. Then the star undergoes 
a catastrophic gravitational collapse that leads to the formation of a black hole, as shown 
schematically in Fig. 6. Here, also the causal light cones, whose generators are the null 
geodesies (192, 193), are depicted. For r < 2M light, and therefore any kind of matter, 
is forced to move inwards towards the singularity. Hence, once the star surface crosses 
the horizon, it becomes disconnected with respect to the rest of the Universe. After this, 
informations on the star cannot any longer be retrieved by observers safely outside the 
horizon. The only information left about the star structure is its gravitational mass M, 
which is also the only parameter in the line element (191) [67]. This is the so called no-hair 
theorem [3]. 

In practice signals from the surface of the star disappear before the formation of the 
horizon. Let one consider the emission of radiation with frequency vq with respect to the 
free falling frame. If the emission occurs at some point with coordinate r (r > 2M), then 
the frequency Voo measured at r> 2M is 

v oo =v {\-—). (194) 

This is a manifestation of the redshift suffered by radiation climbing the gravitational field 
of the star. For r — > 2M, i>oo — » exponentially, namely 

^oc ~ v e~ ku , (195) 

where u is a retarded time defined as 

u — v — 2r* =t-r* (196) 

which is constant along outgoing radial null geodesies (see Eq. (193)): «^coasr^ 2M. 
k is the so-called surface gravity at the black hole horizon. For a Schwarzschild black hole 

AM v ; 

The surface gravity measures the strength of the exponential falloff. It is formally defined 
through the formula [35] 

(X"x*)i#. = ~2fc X M , (198) 
where x M = ( Ji) M ls the timelike Killing vector field whose norm vanishes at the horizon. 
This quantity, as we have seen, plays a central role in quantum effects. 
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Figure 6: Representation of light cones outside a collapsing star (shaded 
region) forming a black hole. Before the horizon formation, and far away 
from the star surface, the light cones are only slightly distorted. After the 
horizon formation, the regions r > 2M and r < 2M (the black hole) are 
causally disconnected: no signal from the r < 2M region can ever reach the 
outer region. 
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The so-called Painleve-Gullstrand coordinates [68] (see the more recent [69]) provide an- 
other regular coordinate system across the future horizon. Moreover they naturally appear 
in the description of acoustic black holes. In ingoing Painleve-Gullstrand coordinates, the 
line element reads 

ds 2 = - (l - dt 2 + ^J^drdt + dr 2 + r 2 (d8 2 + sin 2 Bdtf) . (199) 

The Painleve-Gullstrand coordinates are related to the usual Schwarzschild coordinates 
by the transformation 

j2M/r 

dtFG = dts + l -2M// r ' ( 20 °) 
This is defined in a coordinate patch that covers the asymptotic region (r — > oo) together 
with the region containing the future horizon and the black hole singularity. 
The Painleve-Gullstrand line element exhibits some special features. For instance, the 
constant time spatial slices are flat, with zero curvature. The spacetime curvature of the 
Schwarzschild geometry is encoded into the time-time and time-space components of the 
metric. 

As we have seen, spherically symmetric vacuum (i.e. satisfying = 0) black holes 
are characterized by the gravitational mass M only. However, there also exist charged 
black holes, characterized, in addition, by an Abelian charge Q (which can be either 
electric or magnetic). These are solutions of the full Einstein equations (168) in which 
the source (i.e. T pv ) is the stress energy tensor for the electromagnetic field. Usually the 
latter is expressed in terms of the four- potential 

FfM/ = d»A v - d v A^ , (201) 

and satisfies the curved version of the Maxwell equations 

= , (202) 

whereas the gravitational field is obtained by solving 

Rr* ~ \R9v.v = 8tv(F^F^ - g J^F pa F pa ) . (203) 

These combined Einstein-Maxwell equations admit a unique solution with spherical sym- 
metry, namely the Reissner-Nordstrom [70], with line element 

ds 2 = -(1 - 2M/r + Q 2 /r 2 )dt 2 + (1 - 2M/r + Q 2 /r 2 )~ 1 dr 2 + r 2 dti 2 , (204) 

while the only non-vanishing component of the Maxwell tensor is 

Frt= Q 

4nr 2 

for the electric case (radial electric field), or 



(205) 
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in the magnetic case (radial magnetic field). Eqs. (204, 206) describe the gravitational 
and electromagnetic fields respectively outside a spherical body of charge Q and mass 
M. For M > \Q\, the Reissner-Nordstrom metric Eq. (204), besides the usual curvature 
singularity at r = 0, has two coordinate singularities at 

r± = M ± sjM 2 - Q 2 . (207) 

The singularities at r± coincide when M = \Q\, while no coordinates singularity is present 
for M < \Q\. 

As before, these coordinate singularities are not physical, they are artifacts of the coor- 
dinates system and the spacetime metric can be regularly extended across them. When 
M > \Q\ the solution describes a black hole with horizon located at r+. For r_ < r < r+ 
light and matter are forced to move inwards. However, as the inner horizon r_ is crossed, 
geodesies can avoid the curvature singularity (r = 0) as they bounce back and propagate 
towards an asymptotically flat region, which is different from the one in which they have 
been generated [64]. The bounce is caused by the electromagnetic contribution, namely 
the Q 2 jr 2 term present in the components of the metric, which causes the gravitational 
field to be repulsive at small values of r. 

The horizon at r+ is an infinite redshift surface characterized by a surface gravity 

k rn = . (208) 

The inner horizon r_ is an infinite blueshift surface, which, unlike the event horizon at 
r+, is unstable against perturbations . The infinite blueshift causes perturbations to grow 
and develop a curvature singularity which prevents any continuation across r_. This 
phenomenon is called mass inflation [71]. 

For M = |Q|, the two roots of Eq. (207) coincide and one has a so-called extremal black 
hole with a unique horizon at 

r h = M , (209) 
characterized by vanishing surface gravity. 

Finally, the Reissner-Nordstrom solution for M < \Q\ has no horizons hiding the curvature 
singularity at r = 0. It describes a so-called "naked singularity". The cosmic censorship 
hypothesis forbids the formation of such an object starting from smooth initial conditions. 
However, no explicit proof of this conjecture proposed by Penrose in 1969 [72] is known 
yet. 

D The backreaction in gravitational black holes 

Spacetime curvature has a non trivial effect on quantum fields propagating on it, inducing 
vacuum polarization and particles creation. In the mean, these are described by the 
quantum expectation values (T^v) of the stress energy tensor operator for quantum matter 
and gauge fields. In principle, (T p „) also contain the one-loop gravitons contribution, 
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which is the expectation value of the quantum version of the stress tensor for gravitational 
waves discussed in Appendix B (see Eq. (182)). 

Since gravity couples to the energy momentum tensor, it is clear that the spacetime 
geometry (i.e. the gravitational field) will be affected by the presence of quantum matter 
fields. 

The secular part of the backreaction of quantum fields on the spacetime geometry is 
described by the semiclassical Einstein equations 

~ \R9^ = 87r(T M „) , (210) 

where the r.h.s. should be computed, in the case of our interest, for an arbitrary spherically 
symmetric black hole geometry. The actual black hole metric is obtained by solving Eqs. 
(210) self consistently for g^. Up to now the explicit form of (T^ u (g)) is not known in 
four dimensions; this makes it impossible to solve Eqs. (210). 

Taking into account the vacuum contribution of the quantum matter fields, Schwarzschild 
spacetime is not a solution of the semiclassical equations (210) being the l.h.s. identically 
zero whereas (T i _ lv (g schw )) 7^ (see for instance the analytical approximations developed 
in [53, 73]). 

A black hole dressed by the quantum contribution of the energy momentum tensor can 
be properly described by the metric [74, 75] 

ds 2 = -e 2 ^- v) {l-2m{v,r)/r)dv 2 + 2e 4 '^ v) dvdr + r 2 d{l 2 , (211) 



where the functions ip and m are determined by Eqs. (210), namely 

4nr{T rr ) =4nre i '{T v r ) . (212) 



dr 

To describe an evaporating black hole, the expectation values on the r.h.s. of Eqs. (212) 
are to be taken in the in vacuum state. (A numerical integration of these equations has 
been carried out in [76]. It gives the entire space-time from the collapse till the end 
of the evaporation process. To obtain a closed set of equations, the mean flux (T M , y ) 
has been taken to be that of a two-dimensional field divided by 4ivr 2 . An analytic self- 
consistent calculation of the same flux has been performed in [77] in the adiabatic regime 
dM/dv <JC M 2 . It confirms that the calculation of the Bogoliubov coefficients made in the 
static geometry offers a good approximation as long as the above adiabatic condition is 
satisfied.) 

The apparent horizon of an evaporating black hole is given by 

r h = r(v, r = 2m) = 2M(v) . (213) 

Near this horizon the leading behavior of {T^ v ) represents an ingoing flux of negative 
energy radiation going down the hole [78]. 
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For r ~ rh the evaporating black hole metric (211) can be therefore approximated as an 
ingoing Vaidya metric 

ds 2 = -(1 - 2M(v)/r)dv 2 + 2dvdr + r 2 dQ 2 (214) 
and the relevant backreaction equation reads 

^=4 7 rr 2 (T;> (215) 

giving the rate at which the horizon shrinks in size. By energy conservation, at infinity 
this should equal the luminosity L radiated by the hole. In the case of thermal radiation, 
L ~ AT 4 , where A is roughly the horizon area. 
For a Schwarzschild black hole of mass Mo, this yields 

L ~ AT 4 H ~ M 2 fc 4 ~ 1/M 2 . (216) 

Extrapolating this result in the case of slow variation, one expects the mass loss of the 
system at infinity to be given by 

(217) 

As the mass decreases because of Hawking radiation the process accelerates leading to a 
time scale for the complete evaporation as 

* ~ Mo . (218) 

For black holes with a conserved charge (Reissner-Nordstrom) the evolution appears quite 
different. 

Arguments similar to the one given for a Schwarzschild black hole lead to approximate the 
charged evaporating black hole metric near the horizon with a charged incoming Vaidya 
metric 

ds 2 = -(l-2M(v)/r + Q 2 /r 2 )dv 2 + 2dvdr + r 2 dn, 2 . (219) 
The asymptotic mass loss is, as before, expected to be 

^ ~ —ATh ~ -r\k 4 rn (220) 

where k rn and r+ are given, respectively, in Eqs. (208) and (207). Suppose one starts from 
a black hole with M 3> \Q\. The Hawking temperature increases as the hole looses mass, 
like in the Schwarzschild case. However, when M drops to ^j|<5| one enters a completely 
different regime in which k rn and hence Th decrease as the mass goes down, approaching 
zero when M = \Q\. So the end state of Hawking evaporation for a black hole with a 
conserved U(l) charge is an extremal black hole, which is reached asymptotically in time 
as one can see as follows. For a near extremal black hole (M ~ \Q\) we can write 



Integration gives 



1 ~ MHQ| ' ^ 
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E Quantum stress tensor and vacuum states in 2 
dimensions 



In order to consider backreaction effects due to Hawking radiation in acoustic black holes 
we need to construct the quantum stress {T^) of the phonons (scalar) fields ipi. Starting 
with the action (21), we consider the near-horizon approximation which is equivalent to 
have a free scalar field propagating in the two-dimensional radial part of the spacetime 
metric. Our reduced scalar field is then described by the classical action 



*xy/=g(W>y (223) 
with stress tensor 

r MI , = -^L-^ = a M ^-i(7 M „(V0) 2 . (224) 

The tracelessness property T£ — is a consequence of the conformal invariance of the clas- 
sical action, whereas the conservation V M T M „ = derives from diffeomorphism invariance. 
Unfortunately these two properties cannot hold simultaneously in the quantum theory 
(see for instance [36, 56]). If, as usual, we choose to maintain the latter the covariant 
quantum stress tensor for the scalar fields acquires an anomalous trace which is entirely 
geometrical and state independent 

OP = M , (225) 
where 7? is the Ricci scalar. Choosing a conformal frame for the 2d metric 

ds 2 = -Cdx + dx~ = e 2p dx + dx~ , (226) 

Eq. (225) is translated into 

(T+-) = ~^d + d-p . (227) 

The remaining components (T±±) can be (almost completely) determined by imposing 
the conservation equations 

<MT±±> + 9±(T+_) - 2d±p(T+-} = . (228) 

The result is 

(T±±) = -J^C 1/2 (C- 1/2 ), ±± + t±{T t ). (229) 

127T 

The ambiguity that remains are the two arbitrary functions t±(x ), which encode the 
information about the quantum state in which the expectation values are taken. 

In a generic conformal coordinate system x ± (Eq. (226)) we can define a quantization 
associated to the choice of the positive frequency modes 

(4-7TW) ' e , (A-KW) ' e . (230) 
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Then, if we choose the state | } to be the vacuum state with respect to the modes (230) 
(usually denoted as \x ± )) the functions t± vanish and we get 



{x^T^ix*)^) = -A-( d± pd ±P -dip) . 



(231) 



Moreover, a change of vacuum state 1^*) — ► \x ) produces a change in the expectation 
values of the T±±{x ± ) components of the stress tensor of the form 



where {x ± ,x ± } is the Schwarzian derivative 



± ± _ d 3 Z ± dx ± 3 ( d 2 ^ ,dx^\ 2 
{X ,X } ~ dx± 3 'dx± 2\dx± 2/ dx±) 



(232) 



(233) 



For an acoustic geometry ds 2 



„a" ' dx dx~ , where 
dz 



c\t + 
c(t- 



V 



[ dz 

J c + v 



(234) 
(235) 



ingoing modes are positive frequency with respect to t, whereas near the horizon outgoing 
modes are positive frequency with respect to X~ — —KT 1 e~ KX . So, evaluating the 
Schwarzian derivative, we have t+ = 0, t- = and the stress tensor reads [21] 



(T++) = 
(T__) = 
(T) = 



h c 2 



I2n 4c 4 



It o JZ 

P_ , 3p_ 

" 2p 4p 2 



(c — V ) + V +VV + 



{vv'f 



(T ++ ) + 

h 
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